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ABSTRACT 

We study the properties of luminous stellar "clumps" identified in deep, high-resolution Hubble Space Tele- 
scope NIC2/F160W imaging at L6 /im of six z ^ 2 star-forming galaxies with existing near-infrared integral 
field spectroscopy from SINFONI at the Very Large Telescope. Individual clumps contribute ^ 0.5%- 15% 
of the galaxy-integrated rest-frame w 5000 A emission, with median of k, 2%; the total contribution of clump 
light ranges from 10% - 25%. The median intrinsic clump size and stellar mass are ^ 1 kpc and ^ lO** M©, in 
the ranges for clumps identified in rest-UV or line emission in other studies. The clump sizes and masses in 
the subset of disks are broadly consistent with expectations for clump formation through gravitational instabil- 
ities in gas-rich, turbulent disks given the host galaxies' global properties. By combining the NIC2 data with 
ACS/F814W imaging available for one source, and adaptive-optics assisted SINFONI Ha data for another, we 
infer modest color, M/L, and stellar age variations within each galaxy. In these two objects, sets of clumps 
identified at different wavelengths do not fully overlap; NIC2-identified clumps tend to be redder/older than 
ACS- or Ha-identified clumps without rest-frame optical counterparts. There is evidence for a systematic trend 
of older ages at smaller galactocentric radii among the clumps, consistent with scenarios where inward migra- 
tion of clumps transports material towards the central regions. From constraints on a bulge-like component at 
radii < 1 - 3 kpc, none of the five disks in our sample appears to contain a compact massive stellar core, and 
we do not discern a trend of bulge stellar mass fraction with stellar age of the galaxy. Further observations are 
necessary to probe the build-up of stellar bulges and the role of clumps in this process. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: structure — infrared: galaxies 



1. INTRODUCTION 

Luminous, kiloparsec-sized "clumps" appear to be a ubiq- 
uitous feature of high-redshift star-forming galaxies. Such 
clumps have long been discussed in the context of high resolu- 
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Law et al. 2007). They have been mostly identified in optical 
imaging from the Hubble Space Telescope (HST), probing 
rest-frame UV wavelengths at z > 1. High resolution HST 
near-infrared (near-IR) observations with the NICMOS 
camera, and most recently with the WFC3 instrument, show 
that clumpy structures can also be present in the rest-frame 
optica l ('e.g.,'Toft et al. ' 2007; Dasyra et al. 2008; Kriek ( 
|2009| lElmegreen et al. 200ja; Bu ssmann et al. 
Swinbank et al. 2010a; O verzier et al. 1 12010; Cameron i 
i2010i) . Similar clumps have been identified in rest-frame 
optical line emission from near-IR integral field spect roscopy 
of several z - 2 disks JGenzel et al. I l2006l l200l Hoili) . 
Observations of strongly lensed z ~ 1 - 3 sources in broad- 
band rest-UV/optical light, rest-optical line emission, and 
submillimeter dust continuum and CO line emission with 
source-plane resolution down to ^ 100 pc have confirmed 
characteristic sizes for stellar clumps and star-formin: 
complexes of ^ 300 pc - 1 kpc (e.g. iStarketal.l 120081 
iJones et al.llMToHSwinbank et al.ll2010bl) . 

Various fines of evidence suggest that these clumps are 
not only associated with distinct components in merging 
or interacting systems. They are found in all types of 
galaxies at z > 1 including disks i dentified as such on the 
basis of their structural properties (Elmegreen & ElmegreenI 
2005a; Elme green et al 2005b. 2007,' 2009a} as well 
as their kinematics dOenzel et al.l 120061 l2008l 120111; 
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iBoumaud etaD 120081: iJones etal.l 120101) . High-redshift 
disks are found to have high gas to total baryonic mass 
fractions of ~ 20% - 80% and high local intrinsic gas 
veloc it y dispersions of ~ 20 — 90 k ms~' (e.g., Erb et al. 
2006"; Forster Schreiber et al. ' '2006' '2009; Genzel et a^, 
2008; Crescietal. 2009; Epinatetal. 2009; La w et aL 
2009; Wright etal. 2009; Daddi et al. 2010; Tacco ni et al. 



2010i) . Numerical simulations of gas-rich turbulent disks 



indicate that massive kpc-sized clumps c an form in- situ 
through gr avitational instabilities (e.g., Noguchi 1999; 
Immeli et al. 2004a b; Bournaud, Elmegreen, & Elmegreen 
20071: Elmegreen, Bournaud, & Elmegreen 2008; 
Dekel. Sari. & Ceverino 2009; Agertz, Tevssier, & Moora 
2009 : ICev erino, Dekel, & Bournaud 2010; Aume ret all 
20101: iGenel et al.i i2010). The Toomre length and mass 



derived from the global properties of clumpy disks at high 
redshift correspond roughly to th e clump sizes and masses 



inferred from observa t ions (e.g ., iGenze 



Bournaud et al. I l2008t lElmegr een et akj 



etal. 



l2008l I20TTI 



2009at IJones etal. 



2010h . In these numerical simulations, and according to 



theoretical arguments, the clumps can migrate towards the 
gravitational center as a result of their mutual interactions 
and of dynamical friction against the host disk, and coalesce 
into a young bulge on timescales of 1 Gyr. They may also 
lose part of their mass, contributing to the redistribution of 
mass and angular momentum and to the growth of the disk. 

This "clump-driven" scenario has gained popularity and 
empirical support in recent years, as an important mecha- 
nism for the early formation phases of the spheroidal and 
disk components of present-day massive galaxies. One of the 
major uncertainties in this scenario is the survival of clumps 
against the disruptive effects of tidal torques and of stellar 
feedback, both radiative and mechanical (Murrav et al. 2010l 
Krumholz & Dekell [20T0I: iGenel et al.l [2010: iGenzel et al 



201 ll) . Characterizing in detail the properties of individual 
clumps in distant galaxies, such as their masses, ages, star 
formation histories, and spatial distribution across the host 
galaxy, is crucial to address this issue but remains obviously 
very challenging because of the sensitivity and high spatial 
resolution required. 

In this paper, we study clumps identified in the rest-frame 
optical emission in six massive star-forming z = 2.1-2.5 
galaxies, for which we obtained deep high-resolution HST 
imaging at 1 .6 ^m with the N1CMOS/N1C2 camera and 
F160W filter The targets are among the rest-UV selected 
sources observed wit h SlNFONl as part of the SINS survey 
( Forster Schreibe r et al. 2006, 2009). In a companion paper 
(iForster Schreiber et al. 201 1 , hereafter Paper 1), we derived 
the structural and morphological parameters of the targets 
from the N1C2 imaging, analyzed the results together with 
the Ha line maps and kinematics from SlNFONl, and com- 
pared the morphological properties with those of other sam- 
ples of massive z^2 galaxies imaged with N1C2/F160W but 
selected very differently and probing a wide range in star for- 
mation activity. All but one of our SINS N1C2 targets ex- 
hibit compelling kinematic signatures of disk rotation. The 
global rest-frame optical surface brightness distributions in- 
dicate disk-dominated morphologies as well, with Sersic in- 
dices in the range n « 0.15-2. In contrast, the kinematics of 
the sixth source unambiguously reveal a major merger, the 
main components of which are well resolved in the N1C2 
imaging. Here, we analyze the prominent kpc-scale struc- 
tures seen in all our targets. The combination of SlNFONl 



and N1C2 data sets enables us to determine whether observed 
clumps represent structures consistent with an in-situ forma- 
tion through disk instabilities or, instead, whether they are as- 
sociated with merging components. We focus on the stellar 
properties of the clumps in order to constrain their evolution 
and role in early bulge formation. Two objects are studied 
in more detail, taking advantage of additional high resolution 
HST optical imaging available for one, and SlNFONl obser- 
vations together with adaptive optics (AO) for the other 

The paper is organized as follows. We present the se- 
lection and global properties of the sample, and summarize 
the N1CM0S/N1C2 observations and data reduction in § |2] 
We describe the identification procedure and derive the main 
properties of clumps in the N1C2 images in § |3] In § |4] 
we focus on detailed case studies of the target with exist- 
ing high-resolution HST optical imaging (SSA22a-MD41) 
and the one with AO-assisted Ha SlNFONl observations 
(Q2346-BX482), and further constrain spatial variations in 
colors, stellar mass-to-light ratios, and stellar ages. The im- 
plications of our results for the formation and evolution of 
clumps are discussed in § |5] The paper is summarized in 
§ |6] Throughout, we assume a A-dominated cosmology with 
Hq = 70/170 kms-'Mpc"', n,^ = 0.3, and = 0.7. For this 
cosmology, 1" corresponds to « 8.2 kpc at z = 2.2. Magni- 
tudes are given in the AB photometric system unless other- 
wise specified. 

2. SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

2.1. NIC2 Target Selection and Ancillary Data 

Table [T] lists the six galaxies studied in this pa- 
per The targets were drawn from the initial sam- 
ple of 17 re st-UV selected objects at z ~ 2 f rom the 
SINS survey dForster Schrei ber et al. 1 120061 12009') caiiied 
out with the ne ar-lR inte gral field spectro meter SlNFONl 
(fEisenhaueretal . 2003; B onnet etal. 112004 at the ESO VLT. 
The galaxies were originally part of the large optical spec- 
troscopic survey of z ~ 1.5 -2.5 candida t es se lected by their 
UnGTZ colors described bv ISteidel et al.l (1200 4). Additional 
multi-wavelength data include ground-based near-lR J and 
band imaging, and space-based Spitzer mid-IR photometry at 
3-8 fim with IR AC and at 24 um with MIPS for the major- 
ity of our targets (lErb et al.ll2006t iReddv et al.ll2010l) . For 
one galaxy, MD41, high resolution optical imaging is avail- 
able from observations with the HST Advanced Camera for 
Surveys (ACS) through the F814W filter (hereafter i^u band- 
pass). 

The choice of our N1C2 targets was primarily driven by 
their kinematic nature along with their high signal-to-noise 
(S/N), high quality SlNFONl data mapping the spatial dis- 
tribution and relative gas motions from the Ha line emis- 
sion out to radii > 10 kpc. Five of the targets (BX663, 
MD41, BX389, BX610, and BX482) were explicitly se- 
lected because of their disk-like kinematics, and the sixth one 
(BX528) because it shows in contrast kinematic signatures 
consistent with a counter-rotating binary merger. The galax- 
ies are among the brighter and larger h alf of th e full SINS 
sample in terms of Ha flux and size (see Paper t). The quan- 
titative kinematic classification was performed through appli- 
cation of kinemetry by Shapiro et al. (2008)), an d detailed dy- 
namical mo deling of th e disks is presented by iGenzel et al. I 
( 12008 1) andl CrescietaLl (12009). 

The SlNFONl data sets of our N1C2 targets, among 
the deepest of the SINS survey, are fully described by 
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iForster Schreiber etal.l (|2006[ 12009b . Briefly, aU galax- 
ies were observed in the K band (targeting the Ha and 
[N II] AA6548,6584 emission lines) using the largest pixel 
scale of 125 mas pixel"', five of them in seeing-limited mode 
with FWHM spatial resolution of k, 0"5 (corresponding to 
« 4.1 kpc) and one of them, BX663, with the aid of AO re- 
sulting in a resolution of 0"39 (or 3.2 kpc). Seeing-limited 
//-band data at the 125 mas pixel"' scale were taken for all 
but BX528 to map the H/3 and [O III] AA4959,5007 emis- 
sion lines, also under very good near-IR seeing of « 0"55. 
For BX 482, AO-assisted K band observations were obtained 
using the intermediate pixel scale of 50mas pixel"', with total 
on-source integration time of 6.8 hours and resolution of 0."17 
(1.4 kpc), very similar to that of the NIC2 1.6 /im imaging. 

2.2. Global Properties of the NIC2 Targets 

Our NIC2 sample is relatively bright in the near-IR. The 
five objects with existing K band imaging have a median 
^s,AB = 21.8 mag and span a range of w 1.2 mag in observed /T 
band as well as in rest-frame absolute V band (K^.ab = 21.1 - 
22.3 m ag, and Mv.ab between -23.1 and -21.9 mag; see 
iPaper J) . SpitzerfMlPS 24 fim measurements that have been 
obtained for four of the targets (BX 389, BX 610, BX 528, and 
BX663) indicate that they are among the 3% most luminous 
at rest-frame 8 /xm w avelengths of UV-selected z ~ 2 galaxies 
(iReddv et al.ll20Tol) . 

In terms of global stellar properties, our NIC2 targets probe 
the acti vely star- forming part of the mass ive z^2 gala xy pop- 
ulation (lErb et a l. 2006; Forster Schreib eretZ1l2009l) . Mod- 
eling of their optical to near-IR spectral energy distributions 
(SEDs), including the new photometry obt ained with NIC2, is 
described in Paper I. The results using the iBruzual & CharlotI 
( l2003h evolutionary synthesis code, assumi ng constant star 
formatio n, solar metallic i ty, and adopting the IChabrierl (120031) 
IMF and ICalzetti et al. I (120001) reddening law, are given in 
Table [T] The stellar masses span over an order of mag- 
nitude, from 7.7 x lO' to 1.0 x 10" M© with median w 
5 X 10'" Mq. The median absolute and specific star for- 
mation rates (SFRs) are w 50 Mq yr"' and 0.7 Gyr"' (with 
ranges of 25-185 Moyr"' and 0.6-24 Gyr"'). The objects 
are moderately obscured, with Ay between 0.6 and 1.2 mag, 
or equivalently E(B-V) = 0.15-0.30 mag. The stellar ages 
range from 50 Myr to 2.75 Gyr; BX528, BX663, BX389, 
and BX610 have ages comparable to the age of the Universe 
at their redshift, indicating that they host mature stellar popu- 
lations. The best-fit models to the SEDs of these four objects 
imply rest-frame {U-B)ab ~ 0.75-0.9 mag col ors, close to 
or just about at the separation adopted by, e.g., iKriek et aTl 
(120091) . between "red" and "blue" objects. In contrast, the two 
youngest objects, MD41 and BX482 are significantly bluer 
with rest-frame (//-/?)ab ~ 0.5 mag. 

Two of our targets (BX663 and BX610) are identified 
as candidate AGN on the basis of their obse rved mid-IR 
SEDs from Spitzer IRAC and MIPS photometry (iReddv etal.l 
l2010y) . BX663 also exhibits spectral signatures of Type 
2 AGN in its integrated rest-frame UV and optical spec- 
tra (Shapley et al.l l2004t lErb et al.ll2o5 6'). In our spatially- 
resolved SINFONI data, the central compact emission peak 
of BX 663 is characterized by a higher [N II]/Ha and a broad 
Ha velocity component associated with the AGN, underneath 
the narrower component dominated by star formation. In con- 
trast, BX610 exhibits no sign of an AGN in the rest-frame 
UV and optical, so that its putative AGN is likely to be very 



obscured and will not affect any aspect of our analysis. Al- 
ternatively, its mid-IR emission properties could be due to en- 
hanced PAH emission around rest-frame 8 with respect to 
its far-IR luminosity and in comparison to local star-forming 
galaxies of the same total luminosity, as indicated by recent 
measurements of the full mid- to far-IR SEDs of high-z star- 
forming galaxies (e. g. Nordon et al. 2010; Elbaz et al.ll2010l; 
Muzzinet all 120101; iR. Nordon et al..,2011; ,S. Wuvtset all 



201 Ih . In the context of this paper, we will only consider 



BX 663 explicitly as an AGN source. 

As described in Paper I, our NIC2 targets do not stand out 
compared to the full SINS sample in terms of global stellar 
and dust properties but they do more so in terms of their ob- 
served Ha properties. This is a consequence of our empha- 
sis on selecting objects for the NIC2 follow-up with highest 
S/N and spatially best-resolved Ha emission from the initial 
seeing-limited SINFONI data sets from the SINS survey. The 
median stellar mass and SFR of our NIC2 targets are com- 
parable to those of the SINS Ha sample (w 3 x 10'" Mq 
and w 70 Moyr"'). However, they are among the brighter 
half in terms of integrated Ha fluxes (uncorrected for ex- 
tinction) with median F(Ha) = 1.9 x 10"'^ ergs"'cm"^ (and 
range from 1.1 to 3.1 x 10"'^ ergs"'cm"^), compared to the 
SINS Ha sample median of /^(Ha) = 1 . 1 x 10""' ergs"' cm"^. 
Our NIC2 targets also lie at the high end of the Ha velocity- 
size distribution (Bouche et al. 2007; Forster Schreiber et aH 
l2009l) . Their intrinsic half-light radii ri/2(Ha) ss 4-5 kpc 
are above the SINS sample median of 3.1 kpc, and the cir- 
cular velocities Vd for the five disks are equal to or larger 
than the median for SINS galaxies of « 180 kms"', and 
up to ~ 300 kms"' 0. The ratio of circular velocity to 
intrinsic local veloci ty dispersion of these disks are in the 
range v,; /o-n ?^ 2-6 dGenzel et al.ll2008t ICresci et al.ll2009t 
IForster Schreiber et al. 1 120091) . which is significantly lower 
than for local spiral galaxies (with typical Vc/ctq ^ 10-20; 
e.g., iDib. Bell. & Burkert 2006; Epinat etal.l 120101) . These 
low ratios suggest comparatively larger gas turbulence and ge- 
ometric thickness, and appear to be a characteristic fea ture of 



early disk galaxies at z ^ 1_3 3 (e.g. Forster Schreiber et al. 
'2006; Wright etal.' '2007; Genzel eLalJ |2008; Stark et aL 
2008; .Cresci etal.i ,2009; ,Epinat et al. I [2009; Jones et al. 
12010) . 

2.3. NICM0S/NIC2 Observations and Data Reduction 

The spatial extent and brightness of our targets make them 
particularly well suited for a detailed study of morphologies 
and a reliable assessment of structural parameters from high 
resolution rest-frame optical continuum imaging (iPaper li) . 
Sensitivity and adequate sampling of the instrumental point- 
spread function (PSF) are essential for these purposes, as well 
as for identifying and characterizing small-scale structure, the 
subject of this paper These goals dictated our observational 
strategy and the procedures applied for data reduction, PSF 
characterization, and analysis of the noise properties, which 
are detailed in Paper I and summarized in this subsection. 

The NICMOS observations were carried out between 2007 
April and 2007 September with the NIC2 camera onboard the 
HST and using the F160W filter with mean wavelength of 
1.61 /im (hereafter //i6o bandpass). NIC2 has a pixel scale 
of 0"075, critically sampling the HST PSF at 1.6 /xm. The 

*^ The merger BX528 h as a lower equivalent of 145 kms"'; see 
IForster Schreiber et al. M2009I) for the deiivation. 
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Hifjo bandpass probes rest-frame optical wavelengths around 
5000 A at the median z = 2.23 of our targets, and corre- 
sponds closely to the SDSS g band (which is about 300 A 
bluer in the rest-frame). At these wavelengths, redward of 
the 4000A/Balmer break, the galaxies' light traces the bulk 
of stellar mass more reliably than the rest-frame UV, which is 
significantly more affected by recent star formation and by 
dust extinction. Each target was observed for four orbits, 
amounting to a total on-source integration time of 10240 s 
with good sampling of the PSF. 

The data were reduced using standard routines of the NIC- 
MOS pipeline within the IRAF environment, complemented 
with custom procedures to optimize the removal of bad pix- 
els and cosmic ray hits, and of large scale residuals across 
the images. The images were drizzled onto a final pixel 
scale of 0."05 pixer^ The PSF was determined empirically 
from the two-dimensional profiles of the four brightest, un- 
saturated, and isolated stars present in our NIC2 pointings, 
and has a FWHM of 0."145 based on a Gaussian profile fit. 
The effective noise properties were characterized from the fi- 
nal reduced images, taking into account deviations from pure 
Gaussian noise due to instrumental features as well as corre- 
lated noise and systematics resulting from the data reduction 
procedure. The 3ct limiting surface brightness of the final 
images is /i(i/i6o.AB) = 23.4 mag arcsec"^, and the 3ct lim- 
iting magnitude in a "point-source aperture" with diameter 
J = 1.5 X PSF FWHM = 0."22 (maximizing the S/N of pho- 
tometric measurements in unweighted circular apertures of 
point-like sources) is w 28.1 mag. The total H\a) band mag- 
nitudes measured in a circular aperture with diameter of 3 " 
centered on each target are reported in Table [T] 

Figure [T] presents the final reduced images of the galax- 
ies with sampling of 0"05 pixel"'. The angular resolution 
of FWHM = 0"145 corresponds to a spatial resolution of 
« 1 .2 kpc at the redshifts of our targets. The geometric center 
calculated as the unweighted mean of the x and y coordinates 
of pixels with S/N > 3 is indicated by a cross in Figure[T] For 
the five disks, it corresponds closely to the d ynamical cen- 
ter based on modeling of the Ha kinematics (iGenzel et al.l 
l2008HCresci etal.l2009h . In //jeo band (rest-frame w 5000 A) 
emission, the galaxies exhibit a diffuse and lower-surface- 
brightness component extending out to radii of 0"8- l."2, 
or 6.5- 10 kpc. Rich substructure, including several bright 
compact emission regions, is seen in all sources and con- 
tributes to their irregular appearance. For BX528, two ma- 
jor emission sources to the southeast and northwest are well 
resolved in the NIC2 images (denoted BX 528-SE and NW, 
with projected separation of 8 kpc), with relative H\a) band 
fluxes that suggest a nearly equal-mass major merger (mass 
ratio of '--^ 1.5 : 1), in agreement with the kinemetry analysis 
of the Ha kinematics jShapiro et al. 1200^ . For the large disk, 
BX 389, the small and faint companion 5 kpc in projection to 
the south of, and at nearly the same redshift as the main part 
of the galaxy signals a minor merger (mass ratio of 10 : 1), 
consistent with the lack of substantial disturbances in the Ha 
kinematics. 

3. PROPERTIES OF CLUMPS IDENTIFIED IN THE REST-FRAME 
OPTICAL 

In lPaper fl, we quantified the overall structural properties of 
our NIC2 targets using both parametric and non-parametric 
methods. We compared the rest-frame optical and Ha mor- 
phologies (as well as the rest-frame UV morphology for 



MD41), and also made a differential comparison with other 
samples of massive galaxies at the same redshift but selected 
using different photometric criteria. The effective radius R^, 
Sersic index n, and ratio of minor to major axis b/a derived 
from single-component Sersic model fits to the H\a) band 
maps of our targets are given in Table [T] Now, we turn to the 
prominent and widely distributed clumpy features that consti- 
tute one of the most striking aspects of the rest-frame optical 
morphologies of our SINS NIC2 galaxies. In what follows, 
we use the term "clump" to refer generically to any small- 
scale feature identified with the method described below. This 
method also recovers the merger components of BX 528 as 
clumps, as well as the southern companion of BX 389 and the 
central peak of the AGN-hosting BX 663. 

3.1. Clump Identification and Measurements 

To identify clumps as objectively as possible, and in the 
same way for a ll our galaxies, we used the IRAF task daofind 
jStetsonl Il987h . This task is optimized for identification of 
compact or point-like sources in crowded fields with variable 
underlying background emission, and is thus appropriate for 
our purposes. This task detects objects by searching for local 
density maxima with specified FWHM and peak amplitudes 
above the local background. We used a detection threshold of 
> 3 (7 (with a being the pixel-to-pixel rms noise) and, given 
the observed extent of the bright clumps, we set the char- 
acteristic scales to 1-1.5 x the PSF FWHM. Because the 
clumps are not strictly point sources and have a range of sizes, 
brightnesses, and brightness contrasts against the host galax- 
ies' "background," the number of clumps depends slightly on 
the exact values for the threshold and size scale. Based on ex- 
perimentation, we found that the choices above were the most 
satisfactory according to a visual assessment, and we used the 
same values for all galaxies. A total of 28 clumps are iden- 
tified, with between 2 and 7 per source, and are indicated on 
the NIC2 maps in Figure [T| They are listed in Table |2] along 
with their projected distance t/pioj from the geometric center of 
the galaxies (coinciding with the dynamical center) and other 
properties derived below. We also calculated the deprojected 
galactocentric distance d assuming the clumps are confined to 
a disk of inclination with respect to the sky plane given by the 
best-fit morphological axis ratio of each galaxy (Table [Tji. As 
seen in Figure [T] the two main components of the interacting 
system BX 528 correspond to the neighbouring clumps #1 and 
#2 (BX 528-SE) and to clump #4 (BX528-NW). The central 
peak of BX 663 is identified as clump #2, and the southern 
companion of BX 389 as clump #4. 

Measurements of the clump sizes and brightnesses in our 
NIC2 images are complicated by the bright and variable back- 
ground from the host galaxy, and by the close proximity of 
neighboring clumps often with relatively few resolution ele- 
ments sampling both clumps and interclump regions. More- 
over, large-scale asymmetries in the global light distribution 
of some of the galaxies affect the accuracy of the clump prop- 
ertie s if usin g the residual maps from the Sersic model fits 
(see iPaper t) . The light profiles extracted in circular aper- 
tures of increasing radius centered on the clump positions 
generally show an upturn at approximately mid-distance be- 
tween bright adjacent clumps. For more isolated clumps or 
for those next to other comparatively fainter clumps, the pro- 
files decrease monotonically but show a break where they 
become dominated by the underlying galaxy light or where 
the fainter neighboring clumps start to contribute more sig- 
nificantly. In all those cases, we took the average bright- 
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ness at the upturn or break of the profiles in annuli of width 
dr= 1-1.5 pixels to represent the local background, and com- 
puted the background-subtracted clump brightnesses within 
the apertures of radius rphot just inside these annuli. For the 
other few clumps with no upturn or clear break in their profile , 
we adjusted rphot so as to maximize the aperture size with- 
out overlapping with those of neighboring clumps, with the 
background also measured in narrow annuli just outside of 
the clump photometric apertures. 

The photometric apertures have a diameter dphot in the range 
« 1.4-3.8X the PSF FWHM, on average 2x, and would 
typically miss a significant fraction of the total flux even for 
compact or point-like sources. Since the observed size of the 
clumps is comparable to the spatial resolution and their in- 
trinsic shape is unknown, we applied a correction factor equal 
to the ratio of total flux (enclosed within a radius of 1 ") and 
the flux within rphot based on the PSF curve-of-growth. The 
average correction is a factor of 1.7, with a range from 1.14 
to 2.09. The background measured over rphot < r < rphot + dr 
also includes a contribution from the clump itself, which is, 
for the narrow annuli used, 10% on average (and in the range 
4%- 14%) of the total flux for a point source. This fraction is 
small compared to both the aperture correction and the back- 
ground estimate, so we neglected this additional correction. 
We determined clump sizes by taking the direct FWHM of 
the background-subtracted light profiles (i.e., twice the ra- 
dius at which the profile is half its maximum value). Nine 
of the 28 clumps are unresolved, with direct FWHM smaller 
than the PSF FWHM; we adopted the PSF FWHM as an up- 
per limit on their size. For the others, the intrinsic FWHM 
size is simply calculated by subtracting in quadrature the ob- 
sei-ved clump and PSF FWHMs. The FWHM sizes, Hm) band 
magnitudes, and fractions of the total galaxy light contributed 
by the clumps (denoted FWHM"^', Hf^^, and T^^) as well as 
the photometric aperture radii and corrections are given in Ta- 
ble |2] 

To gauge possible uncertainties on the clump light contri- 
butions, we considered two other measurement methods. In 
the first, we used the peak flux above the local background 
as returned by daofind when identifying the clumps and com- 
puted the total flux assuming Gaussian unresolved sources. 
Since the clumps are mostly (s^ 2/3) consistent with being 
extended — albeit compact — sources, this method is ex- 
pected to provide lower limits. The inferred fluxes are on 
average (and median) a factor of « 3 lower than the adopted 
estimates described just above, and with smaller differences 
typically found for the unresolved clumps (with FWHM"^' < 
PSF FWHM). Alternatively, we estimated the integrated light 
within apertures of radius rphot but without background sub- 
traction and apertur e correction, in a manner analogous to 
that used by Elmegr een et al. I (l2009al) . In view of the signif- 
icant surface brightness of the underlying light in our galax- 
ies, these estimates represent upper limits and are typically a 
factor PS 2.5 higher than the adopted values. The correspond- 
ing fractional light contributions obtained with these two ap- 
proaches ( J^p5 and F^l^, respectively) are also given in Table|2] 
for comparison with the adopted measurements. We thus in- 
fer that the adopted clump light contributions are uncertain to 
a factor of ss 3. 

We emphasize that the contrast between the clumps and 
the surrounding, or interclump regions in the i/ieo band im- 
ages of our NIC2 targets is relatively small. In the photo- 
metric apertures employed, the total contribution from the lo- 



cal background light is typically 3 - 4 times higher than 
the background-subtracted clump fluxes, making the mea- 
surements of the clump brightnesses rather sensitive to the 
background estimates as seen above. This is very different 
from the "clump-cluster," "chain," and spiral galaxies from 
flie H ubble Ultra Deep Field (H-UDF ) studied in the opti- 
cal bvlElmegreen & Elmegreen! (I2005al) andlElmegree n et al. I 
(I2009al) . in which the clumps identified in ACS F775W imag- 
ing (;775 band) are typically brighter by factors of ^ 2-4 than 
the surrounding emission, with the contrast increasing fur- 
ther a t shorter wavelengths. While Elmegreen & Elmegreen] 
(l2005a) accounted f or the galaxy's ba ckground in measuring 
clump magnitudes, lEknegreen et al. I ( 12009 a) neglected this 
contribution since it is generally small in the optical for their 
sample of H-UDF galaxies. 

We computed the approximate observed (i.e., uncor- 
rected for dust extinction) rest-frame SDSS g band lumi- 
nosities of the clumps from the background-subtracted and 
aperture-corrected //i6o band magnitudes, adopting, to ex- 
press the results in solar units, an absolute Mg = 5.07 mag 
dBlanton eT al. 2001)E3- We further estimated the clump stel- 
lar masses by multiplying the My, derived from the SED mod- 
eling of each galaxy with F'^'^. This scaling corresponds to 
assuming, crudely, a constant ratio across each source be- 
tween stellar mass and observed H\a) band light or, equiva- 
lently, rest-frame g band light. For MD 41, the M^./L'^'*' ratio 
map obtained from the /sw-Z/ieo color map constructed from 
its NIC2 and ACS images indicates a fairly small scatter of 
~ 0.2 dex in observed log(Mi/L^''^') across the galaxy (see 
iPaper I[) . For the clumps themselves, the more detailed anal- 
ysis presented in §|4]for MD 41 and BX482 suggests that the 
assumption of constant M/L ratio is on average valid to within 
a factor of ^ 2 - 3. The clump luminosities and stellar masses, 
^rest.ci ^jjj ^ci^ ijgjgjj Table El 

3.2. Clump Properties for the NIC2 Targets 

Figure |2] plots the //160 band magnitude, fractional light 
contribution, and inferred stellar mass of the clumps as a 
function of observed size. The distribution of intrinsic sizes 
and stellar masses of the resolved clumps is shown in Fig- 
ure [3] The observed FWHM sizes range from unresolved up 
to « 2 kpc. Excluding the nine unresolved clumps, the in- 
trinsic FWHM" sizes span « 0.1-1.5 kpc, with mean and 
median of w 0.9 kpc, or ss 0"1 1. The clump magnitudes vary 
between ^^J^goAB ~ ^4.3 and 27.9 mag, with mean and me- 
dian of w 26.5 mag. On average, the clumps have « 4% 
(median w 2%), ranging from approximately 0.5% to 16%. 
With the constant M/L ratio assumption for each source, these 
fractions imply Mf ^ 9 x 10^ to 9 x 10'^ M©, and typically 
^10' Mq. Considering only the resolved clumps, the mean 
and median Mj' ^ 2 x 10^ Mq. Among the resolved clumps, 
we find no significant correlation between clump magnitude 
or fractional light contribution and size, and the scatter of the 
data is fairly large (« 0.9 mag in //leo, and about equal to 
the mean in J^*^'). There is also no obvious trend of clump 

The NIC2 //uo bandpass probes a wavelength range between the rest- 
frame g and V bands for our targets, with closer effective wavelength and 
larger overlap with the g bandpass. For the best-fit SED parameters of the 
sources, the color term between observed g and V band magnitudes, or lumi- 
nosities, varies from 3% to 20% depending on the object, with median/mean 
of 10%; the difference between the actual rest-frame luminosity probed 
through the //|60 bandpass and the g band would be smaller 
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size, magnitude, fractional light contribution, or stellar mass 
as a function of galactocentric distance (c/proj or d) nor of 
distance normalized to the R^ff of the galaxies. In Figure |2] 
panel a, dotted curves indicate constant apparent central sur- 
face brightnesses, i.e., within a circular aperture of diameter 
equal to the observed FWHM'^'. The clumps cover about 2 
magnitudes in central surface brightness, with a standard de- 
viation of « 0.7 mag arcsec"^ around the median (and mean) 
value of /ii6o.AB ~ 23.7 mag arcsec"^. This value corresponds 
roughly to the 3 a surface brightness limit of our NIC2 im- 
ages. 

The "clumps" associated with the merger components of 
BX528, the southern companion of BX389, and the central 
core of BX 663 (plotted with different symbols in Figures |2] 
and|3]l do not differentiate in terms of size but are part of the 
brighter half and have J^"^' above the median value. Excluding 
them from our clump sample does not, however, significantly 
affect the ranges and the mean or median properties inferred 
above nor the absence of correlation between clump proper- 
ties. 

Despite being a striking feature in the appearance of all our 
targets, our results imply that the individual clumps contribute 
typically only a few percent of the total light. Even with the 
substantial uncertainties of a factor of « 3, each clump still 
represents a fairly small contribution. Summed up over all 
clumps identified in a given galaxy, the total contributions 
are typically ^ 15% per object, ranging from roughly 10% 
to 25% (at least ~ 5% from the lower limits limits obtained 
in the alternative "point-source" estimates made above, and at 
most ^ 50% from the upper limits derived without subtract- 
ing the host background in the photometric apertures). Thus, 
for our targets, the bulk of the rest-frame « 5000 A light ap- 
pears to be dominated by an underlying component from a 
more smoothly distributed stellar population. Alternatively, 
fainter, unresolved clumps and clusters may also contribute to 
this extended emission component. 

We searched for correlations between clump properties and 
the global properties derived from the SED modeling of the 
host source. There is no significant correlation between clump 
size, magnitude, or fractional light contribution with galaxy 
stellar mass, age, absolute and specific SFR, nor between 
clump size and fractional light contribution and extinction. 
The only possible trend is of fainter clump //i6o-band magni- 
tude with higher Ay, with Spearman's rank correlation coeffi- 
cient p = 0.53 and 2.8 a significance level The trend weak- 
ens when excluding the clumps associated with the merging 
components of BX528, the southern companion of BX389, 
and the central peak of the AGN source BX663 (p = 0.37, 
1.7 cr). There is no apparent trend with morphological pa- 
rameters although it is noteworthy that BX610 and BX482, 
the two disks with qualitatively most prominent clumps, have 
higher Gini coefficients than MD41 and BX389, the two 
disks where the clumps appear visually less regular in shape. 
MD 41 and BX 389 are the most edge-on disks based on their 
axis ratio and have the largest inferred Ay, so that higher and 
possibly patchy extinction may affect the appearance of these 
galaxies more than for the others. Our sample is small, how- 
ev er, a nd covers a limited range in host galaxy parameters 
(§ I2.2I 1. so that it would be difficult to discern all but the 
strongest possible trends. 

The significance quoted corresponds to tlie number of standai'd devi- 
ations by which the sum-squared difference of ranks deviates from the ex- 
pected value under the null hypothesis of statistical independence (p = 0). 



3.3. Constraints on Emission Line Contributions 

At the redshift of our targets, the NIC2 F160W bandpass in- 
cludes the [O III] AA4959,5007 doublet and H/3. It is thus im- 
portant to assess whether the observed //leo band morpholo- 
gies are biased by nebular line emission. For all sources ex- 
cept BX 528, SINFONI seeing-limited observations of these 
lines w ith the H band grating were taken as part of the SINS 
survey (Forster Schreiber et al. 2009), and are analyzed in de- 
tail by P. Buschkamp et al. (2011, in preparation). The [OIII] 
lines are detected in all galaxies except BX663, and H/3 is 
additionally detected in MD41 and BX610. Some of the 
non-detections can be attributed to the galaxies' lines being 
redshifted to the wavelength of bright telluric lines, in par- 
ticular H/3 for BX663 and BX482. Based on the integrated 
line flux es (or the 3 cr upper limits thereof), we concluded in 
iPaper j that the broad-band Hia) emission of our NIC2 tar- 
gets is dominated by the continuum, with only a modest to 
very small fraction from the [O III] and H/3 lines. The largest 
line contribution is inferred for BX 389, where the [O III] lines 
make up 24% of the integrated //iso band flux and, with the 
3(7 upper limit on H/3, the total line contribution is estimated 
at < 29%. BX610 has the smallest contributions, amounting 
to w 6% for all three lines combined. Our examination of the 
sources on a pixel-to-pixel basis further indicates no impor- 
tant spatial variations in the emission line contributions. 

In the context of this paper, it is worth considering more 
specifically whether the clumps identified in § 13. II could be 
due to nebular line emission rather than stellar continuum 
emission, since their contribution also represents a modest 
fraction of the total //160 band light. A detailed assessment 
would require AO-assisted SINFONI H band observations, 
but the existing seeing-limited data provide useful indications. 
For BX610 and BX482, the total fraction of light from the 
identified clumps (16% and 25%, respectively; Table|2j is at 
least twice a s high a s the total emission line contribution (6% 
and < 12%: lPaper J) , implying a major contribution from the 
continuum in these two sources. For BX663, MD41, and 
BX 389, the line contributions (or 3 cr upper limits thereof) are 
comparable to or higher than the sum of the clump contribu- 
tions so that it is not possible to set the same constraints with 
the current data. For BX482, the spatial distributions in Ha 
and i/i6o band follow each other fairly well, albeit with some 
small differences discussed in the next section. Our analysis 
in §|4]implies that the Ha to Hi^) band flux ratio at the loca- 
tion of the clumps identified in the NIC2 image is a factor of 

< 2 the source-integrated ratio. If we assume approximately 
constant observed emission line ratios across the galaxy, this 
would imply [O III] and H/3 contributions to the //160 band 
emission of at most twice that of the global contribution, i.e., 

< 25% (Paper I). Again, this means the clumps need to have 
a major component of continuum light in Hi^) band. From 
the estimates above, and at least for BX610 and BX482, 
the clumps in the NIC2 images do not appear predominantly 
caused by line emission. More generally, clumps identi- 
fied in high resolution optical and near-IR broad-band imag- 
ing of other samples (e.g., lElmegreen & Elm egreen 2005^ 
lElmegreen et al. 2009a) are ubiquitous in z > 1 galaxies, over 
a range of redshift that precludes an origin entirely from emis- 
sion lines in all cases. 

4. DETAILED CASE STUDIES: MD 41 AND BX 482 

The sizes, luminosities, and stellar masses derived in the 
previous section reflect the current state of the clumps but 
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provide little information on their evolution. In addition, the 
stellar masses were based on the crude assumption of constant 
M /L ratio across the galaxies. In the following, we examine 
two of our targets with existing high spatial resolution ancil- 
lary data to assess possible variations in M/L ratios and search 
for differences in age-sensitive properties among clumps, and 
between clump and interclump regions. We combine our 
NIC2 maps with the available ACS imaging for MD 41 to look 
at the !8i4-//i60 colors, and with the AO-assisted SINFONI 
Ha map of BX 482 to determine the Ha equivalent widths. As 
we will see, the /su-^^ieo colors provide notably tighter con- 
straints on the M/L ratios, irrespective of dust extinction, stel- 
lar age, and star formation history. On the other hand, the Ha 
equivalent widths are insensitive to extinction if the line and 
continuum emission originates from the same regions, and are 
thus most useful to constrain variations in stellar ages. 

4.1. MD41: Constraints from NIC2 and ACS Observations 
4.1.1. Rest-frame Optical- and UV-Identified Clumps 

The ACS ;8i4 band map of MD41 is shown in Figure |4] 
alongside the NIC2 //leo band map and a color-composite im- 
age. For the analysis, the ACS map has been convolved to 
the spatial resolution of the NI C2 data and registered onto 
the same astrometric frame (see iPaperll for details). The ap- 
pearance of MD41 is broadly similar in igu and //leo band. 
The measurements in iPaper II of the structural and morpho- 
logical parameters yield comparable effective radius, Sersic 
n index, axis ratio, as well as Gini G, Multiplicity ^I*, and 
M20 coefficients for both images. The /gi4-//i6o colors across 
MD41 are also fairly uniform, with median of 0.95 mag and 
standard deviation of 0.37 mag over the regions where pix- 
els have S/N > 3 in both bands. These results indicate that 
there is no strong A:-correction between the rest-frame UV and 
optical morphologies for MD41, consistent with the general 
findings for rest-UV selected samples at z 2, and in con- 
trast to those for samples selected at longer wavelengths with 
typically redder rest-frame optical colors (iPeter etal. 1 120071 : 
iLaw et al.ll2007l:lToft et al.llIOOTtlCameron et al.ll20I0 l). 

However, it is also evident from Figure HI that the bright- 
est peaks do not all coincide spatially between the two maps. 
Therefore, we applied the procedure described in § B.ll to de- 
tect clumps in the ACS image matched to the PSF of the NIC2 
data, using the same detection parameters in the IRAF task 
daophot. Four clumps are thus identified in /gi4 band emis- 
sion, one of which corresponds to NIC2 clump #4 whereas 
the three others have no counterpart among the NIC2 clumps 
(hereafter referred to as 'ACS-#1," 'ACS-#2," and "ACS- 
#3")- All of the clumps are indicated on the maps of Figure|4] 
We computed the clump photometry with background sub- 
traction and aperture correction, as well as without for com- 
parison purposes (hereafter simply "background-subtracted" 
and "raw" photometry). We measured the /814 band bright- 
ness of the NIC2 clumps in the same apertures used to de- 
rive their i/ieo band brightnesses. For the three additional 
ACS clumps, we defined appropriate apertures as explained in 
§ l3.1l but based on the /814 band light profiles. The ACS-only 
clumps have no local excess emission above our detection 
threshold in the NIC2 map, so we adopted 3 a upper limits 
on their background-subtracted Hia) band magnitudes based 
on the noise properties of the image. For the NIC2-identified 
clumps, the background-subtracted /814 band magnitudes are 
formally above the 3 cr noise level in the photometric apertures 
except for clump #1, for which we used a 3 tr upper limit. 



4.1.2. Variations in Colors of the Clumps 

The photometric measurements are reported in Table |3]and 
plotted in Figure|5] Panel (a) shows the //if,o versus /gu band 
magnitudes for individual pixels with S/N > 3 in each band. 
Panel (b) shows the magnitudes estimated for the clumps, 
with and without subtraction of the galaxy's background. For 
reference, the total magnitudes of MD41 (in a circular aper- 
ture of 3" diameter) and those corresponding to the sum of 
the fluxes of all pixels with S/N > 3 in both bands are shown 
as wellQ In both panels of Figure |5] the diagonal lines cor- 
respond to constant /gu-Z/ieo colors. A vector shows the 
significant impact of extinction in the rest-frame of the galaxy 
on the observed colors. 

AU the NIC2 and ACS clumps in MD41 contribute small 
fractions of ~ l%-5% to the galaxy-integrated fluxes in 
both //160 and /814 band emission (rest-frame w 5000 A and 
« 2600 A), respectively. Overall, the clump and interclump 
regions in MD 41 span a comparable range in their colors, re- 
flecting the uniformity across the galaxy noted above. The 
median color of the interclump regions is 1.0 mag, nearly 
equal to the median over all S/N > 3 pixels (0.95 mag). The 
colors of the clumps are mostly within « 0.5 mag of these val- 
ues, for both the raw and background-subtracted estimates (or 
limits). Interestingly, there appears to be a systematic differ- 
ence in colors between the clumps identified in the different 
bands: the NIC2 clumps have a median color roughly 0.6 mag 
redder than the ACS clumps. The two most distinct clumps 
are NIC2-#1 and ACS-#3, which have colors up to 2 mag 
redder and bluer, respectively, than the ensemble of the other 
clumps and the interclump regions. 

4. 1.3. Implications for the Stellar Properties of the Clumps 

With the /814-//160 colors of the clumps in M D41, we can 
assess the reliability of the assumption made in § 13. li on their 
M/L ratio and the consequences on the clump stellar masses. 
To interpret the colors, we used the me an /gw-Z/ieo versus 
M-^/U^^^ relationship derived in iPaper II from the evolution- 
ary synthesis models of Bruzu al & Chariot (l2003h for a wide 
range of star formation histories (SFHs), stellar ages, and dust 
extinction. The observed /gu-Z/ieo colors are strongly de- 
generate in age, extinction, and SFH. On the other hand, the 
various model curves occupy a fairly well defined locus in 
the ii\A-H\a) versus Mi^/U^^^ parameter space, and the stan- 
dard deviation of the models suggest the relation is accurate 
to 0. 15 dex in log(M^/L'^''^') over the range of colors observed 
across MD 41, with little sensitivity to metallicity between so- 
lar and 1/5 solar (appropriate for our NIC2 targets includ- 
ing MD41: lP. Buschk amipet al.l l2011. in preparation!) . Stel- 
lar masses were then estimated by multiplying the M/L ratios 
with the observed luminosities from the //160 band data. The 
color-based M^/U^^^ and estimates for all clumps identi- 
fied in MD 41 are given in Table [3] 

The color distribution of individual pixels across MD41 
implies a median log(M*/L^'''*' [M© L© g]) « -0.8 over the in- 
terclump regions (-0.85 over all S/N > 3 pixels). For most 
clumps, the ratios (or limits thereof) derived from the colors 
agree within a factor of ~ 2 - 3 with the median over the in- 
terclump regions. The most notable differences are for the 

" The small difference between these two "integrated" measurements is 
due to the different areas sampled, with the total aperture including the two 
sources ^ 1" southeast and northwest of MD41 seen in the ACS image but 
undetected in the NIC2 image. 
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reddest and bluest clumps, NIC2-#1 and ACS-#3. The ra- 
tio between the Mi,/U^^^ values (or limits thereof) for these 
two clumps ranges from « 8-65, depending on the treatment 
of the galaxy's background. The systematic offset in colors 
between the NIC2- and ACS -identified clumps implies over- 
all higher Mi,/U^^^ ratios and larger M^, for the former set of 
clumps by factors of '--^ 2 and ~ 4 - 6, respectively (or more 
considering the limits). For the NIC2-identified clumps #2, 
#3, and #4, the background-subtracted color-based masses in- 
ferr ed he re are k, 1.5-2.5 times larger than those calculated 
in § l3.1l with the assumption that clumps have the same M /L 
ratio as the entire galaxy; for clump #1, it is almost a factor 
of eight higher The comparison of the trends based on the 
raw versus background-subtracted photometry in MD41 in- 
dicates that, while still fairly modest between most clumps, 
the variations in colors, hence M/L ratios tend to be accen- 
tuated (along with a general shift towards redder colors and 
higher M/L ratios) when accounting for the galaxy's back- 
ground light contribution. 

The observed distributions of colors in MD 41 could be in- 
terpreted as indicative of overall modest variations in extinc- 
tion and ages between clumps and across the galaxy, with 
the apparent exceptions of NIC2-#1 and ACS-#3. The dif- 
ference in median /gw-Z^ieo color between the NIC2- and 
ACS -identified clumps by 0.6 mag could be accounted for by 
a difference in Ay of about 0.7 mag. Relative to the inter- 
clump regions, the redder color of NIC2-#1 and bluer color 
of ACS-#3 by up to « 2 mag would require differences in Ay 
of up to about 2.5 mag; given the best-fit Ay = 1.2 mag for 
the integrated SED of MD41 (Table [U, it appears plausible 
that not only extinction but also age variations contribute to 
the observed clump colors, at least for these two clumps. If 
interpreted solely in terms of stellar age, assuming constant 
star formation (CSF) and a uniform Ay = 1.2 mag, the NIC2- 
only clumps would be on average ~ 5 - 1 times older than the 
ACS-only clumps. Shorter timescales for the star formation 
activity would tend to reduce (but not eliminate) the overall 
relative age difference between NIC2- and ACS-only clumps, 
if extinction variations are neglected. 

The results for MD41 suggest that, with the possible ex- 
ception of the reddest and bluest clumps, the M/L and stellar 
masses are reliable to a factor of a few. The NIC2-identified 
clumps tend to be redder and more massive than the ACS- 
identified clumps. We find indications that not only extinction 
effects but also variations in stellar age are likely responsible 
for the full range in colors of the clumps. The main uncer- 
tainties affecting our analysis are associated with the limited 
photometry available for MD41 and with the difficulties in 
evaluating the clump light in the presence of a bright back- 
ground from the host galaxy. 



4.2. BX482: Constraints from NIC2 and SINFONI 
Observations 

4.2.1. Rest-frame Optical- and Ha-Identified Clumps 

For BX 482, we can set more robust constraints on relative 
stellar ages between the clumps and the interclump regions 
with the help of the SINFONI AO-assisted Ha map available 
for this galaxy. Ideally, the Ha line flux combined with the 
continuum emission measured in the near vicinity of Ha in 
wavelength would provide the Ha equivalent width, which 
varies strongly with the age of a stellar population but is insen- 



sitive to dust extinction^ However, the continuum map from 
our SINFONI data of BX 482 has insufficient S/N for our pur- 
poses, so we used the H\(,q map as a proxy. In what follows 
below, the Ha equivalent widths are calculated in the rest- 
frame (denoted W'^^'(Ha)) assuming a flat f,y continuum flux 
per unit frequency between rest-frame 5000 A and 6563 A and 
using the wavelength of Ha for the conversion to continuum 
flux per unit wavelength, f\ . 

We analyzed the NIC2 //leo and SINFONI Ha data of 
BX482 in the same fashion as the NIC2 and ACS data of 
MD 41, matching the NIC2 PSF to that of the SINFONI data, 
which has a FWHM = 0"17. The NIC2 and Ha images were 
aligned relying primarily on the outer isophotes and on the 
brightest clump visible in both maps (NIC2-#1). We refer the 
reader to Paper I for a full description of these procedures. 
The NIC2 H\(,q band and SINFONI Ha maps are shown in 
Figure|6]along with a color-composite image created from the 
PSF-matched maps. Both the lower surface brightness com- 
ponent and the brighter substructure are overall very compa- 
rable between the //i6o band and Ha emission. This similarity 
is reflected in the nearly equal values of G, ^I', and M20 coeffi- 
cients, and the structural parameters /?eft, n, and b/a from both 
images agree well within the 1 a uncertainties (see PaperT). 

However, some differences do exist in the details of the 
spatial distributions. Specifically, clump #4 identified in 
the NIC2 data does not have an obvious counterpart in Ha, 
whereas there are peaks in Ha, significantly offset to the north 
of clump #4 and on the northwestern side, that are not appar- 
ent in the NIC2 image. We thus ran the clump detection pro- 
cedure on the Ha map, finding matches for all NIC2 clumps 
except #4, and two additional clumps in line emission, here- 
after denoted "Ha-#1" and "Ha-#2." The clumps are indi- 
cated on the maps of Figure |6] We measured the H\(,q band 
and Ha fluxes of each clump, within the same aperture for 
consistent photometry, and both with and without subtraction 
of the galaxy's background. An aperture correction was ap- 
plied for the background-subtracted measurements. The same 
apertures as in § 13 . 1 I were used for the NIC2-identified clumps 
(the H\a) band magnitudes measured from the original and 
PSF-matched NIC2 maps are not significantly different). For 
the Ha-only clumps, the apertures were chosen based the Ha 
light profiles. For clumps undetected in one or the other map, 
we adopted 3 a upper limits on their background-subtracted 
fluxes, derived from the noise properties of the images. The 
flux measurements are given in Table |4] 

4.2.2. Variations in Ha Equivalent Width of the Clumps 

Figure |7] shows in panel (a) the distribution of H\(,q band 
flux densities versus Ha line fluxes for individual pixels with 
S/N > 3 in each map. In panel (b), the results for the 
clumps (with and without background subtraction) are plot- 
ted along with the galaxy-integrated measurements and the 
fluxes summed over S/N > 3 pixels for comparisonPl Di- 

This statement is true unless different amounts of dust obscure tlie 
H II regions compared to the bulk of the stellar population, or a significant 
fraction of Lyman continuum photons is absorbed by dust within the H II 
regions or escapes from the H II regions. In the analysis presented here, 
we neglected these effects. There is eviden ce for such "differential attenu- 
ation" on global galactic scales locally (e.g.JCalzetti et al. 2000) and a t high 
z (e.g ., [Forster Schreiber et al. 2009; Mancini et al. 20 1 1 : Yoshika wa et al. I 
1201 Vti but if it is similar at different locations across galaxies, the relative age 
estimates would not be significantly affected. 

^' The small differences between the latter two is due to fewer pixels 
meeting the S/N criterion on the fainter northwest side than on the southeast 
side within the radius of \'.'5 of the total photometric aperture. 
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agonal lines are overplotted in both panels, corresponding to 
constant values of W^'^^'(Hq;) and labeled with equivalent stel- 
lar ages for CSF models as described below. Extinction does 
not have a significant effect on the W"^'^'(Ha) due to the prox- 
imity in wavelength of Ha and the rest-frame w 5000 A re- 
gion probed by the //leo bandpass; the effects of extinction run 
nearly parallel to the lines of constant Ha equivalent width in 
Figure |7] 

Most clumps in BX482 contribute ^ l%-5% of the to- 
tal emission in both //igo band and Ha, and the clumps and 
interclump regions span a comparable range of Ha equiva- 
lent widths. While there are some differences for individual 
clumps between the background-subtracted and raw flux es- 
timates, the median properties are almost the same. In the 
ensemble, the clumps have higher W^''^'(Ha) than the inter- 
clump regions, with median values of 210 A and 150 A, re- 
spectively. The two Ha-identified clumps have the highest 
background-subtracted equivalent widths. Two clumps are of 
particular interest in BX482: NIC2 clump #1 and #4. Clump 
#1 is the largest, most luminous, and one of the most massive 
of all clumps identified in ourNIC2 targets (see § 13.2) . It also 
has the largest contribution to the total host galaxy's emis- 
sion, ^ 15% in Hia) band as well as in Ha, roughly 5 - 10 x 
higher than the average clump in our NIC2 targets. How- 
ever, clump #1 does not differentiate in terms of its W^'''^'(Ha), 
which is equal to the median value of all clumps in BX482. 
In contrast, clump #4 has about the same size and fractional 
light contribution as the median values for the clumps in our 
NIC2 targets, but it is one of the least massive clump. It is 
further distinct among the BX482 clumps, with the lowest 
measured W^''^'(Ha) ^ 90- 120 A, depending on whether the 
galaxy background is subtracted or not. 

4.2.3. Variations in Ages of the Clumps and Implications for the 
M/L Ratios 

To interpret the Ha equivalent widths in terms of 
stellar ages, we compu ted model predictions using the 
iBruzual & Chariot* ('2003') evolutionary synthesis code with 
a Chabrier (2003) IMF and solar metaUicity. The Ha 
flux was calculated from the H ionizing rates of the mod- 
els, applying the recombinatio n coefficients for case B 
from iHummer & StorevI (1 19871) for an electron tempera- 
ture Te = 10^ K and dens ity «e = 10^ cm""* (see, e.g., 
iForster Schreiber et al. 11200^ . and the //i6o-band flux density 
was approximated from the g-band luminosity synthesized 
from the model spectra. The effects of stellar photospheric 
absorption on our Ha flux measurements are negligible, and 
so were not accounted for in the model predictions. The vari- 
ations of VK'^'^'*'(Ha) with stellar age depend on the SFH. The 
lines of constant W^''*"(Ha) plotted in Figure|2]are labeled with 
the corresponding ages for a CSF model. As an example of 
alternative SFHs, for exponentially declining SFRs with e- 
folding timescales r, the model curves for the same ages in 
Figure |7] shift downwards along the vertical axis and become 
more widely separated for ages > t. 

For CSF, the W^'''^'(Ha) measurements imply a median stel- 
lar age of 55 Myr for the clumps and 140 Myr for the in- 
terclump pixels. For the large, luminous, and massive NIC2 
clump #1, its W"'^'(Ha) also implies an age of 55 Myr, and 
so it does not appear to stand out in its evolutionary stage 
compared to most other clumps. The smaller, low-mass 
NIC2 clump #4 with lowest W"^'^'(Ha) is the oldest clump in 
BX482, with inferred age of ~ 240-570 Myr depending on 



the treatment of the background galaxy light. The impact of 
various SFHs, with SFR declining or increasing with time, on 
the derived ages of the clumps and interclump regions is illus- 
trated in Appendix IaI In all cases, however, the relative ages 
are qualitatively unchanged as long as the different regions 
have similar SFHs. It is conceivable that the duration of the 
star formation activity is shorter for the localized kpc-sized 
clumps compared to the bulk of the stellar population across 
the galaxy. If so, most clumps would still be younger than the 
interclump regions unless they follow exponentially increas- 
ing SFHs with timescales t < 50 Myr (see AppendixlAll. 

We can verify indirectly the validity of our hypothesis in 
§ l3.1l of constant observed M/L for BX482, on the basis of 
the relative ages and assuming that extinction plays no role 
(i.e., that the relative dust obscuration is always the same 
across the galaxy and bet ween the clumps). Using the same 
IBruzual & Chariot! (12003") models, and for a given SFH, the 
ages inferred for the clumps and the interclump regions im- 
ply comparable M/L ratios: the scatter among clumps is 
< 0.3 dex, and the median log(M*/L^'''*') values of the clumps 
are all within 0.2 dex of each other for different SFHs. 
These results for BX 482 provide another example where the 
M/L of the clumps appears to be roughly constant and compa- 
rable to the ratio across the entire galaxy, as we inferred from 
the analysis of the /sw-i/ieo colors for MD41. Obviously, 
this conclusion should not be generalized to all galaxies, but 
it supports the plausibility of such an assumption in the ab- 
sence of better empirical constraints. 

The analysis of the extinction-insensitive Ha equivalent 
widths in BX482 presented in this section revealed signifi- 
cant variations in the relative stellar ages among the clumps. 
If similar SFHs are applicable, the data provide indications 
that all clumps have stellar populations of comparable or 
younger ages than the interclump regions, with the excep- 
tion of one NIC2-identified clump without Ha counterpart. 
In addition, there is no evidence for variations by more than 
a factor of a few in observed M/L ratios between clumps and 
interclump regions. The main uncertainty is in the SFH ap- 
propriate for different regions in BX482, which cannot be 
constrained with the data currently available. In future stud- 
ies, detailed spatially-resolved constraints on the SFHs across 
galaxies will help in this respect. More accurate knowledge of 
absolute ages will also be important to constrain clump life- 
times, a crucial issue in the context of clump-driven scenar- 
ios for the evolution of tu rbulent, gas-rich dis ks at high red- 
shift (see the discussion by lGenzel et al.llToi IL and references 
therein). 

5. DISCUSSION 

Our actively star-forming, rest-UV selected NIC2 targets 
exhibit clumps in their rest-frame optical emission, just as 
clumps are observed in the rest-frame UV emission of many 
z > 1 star-forming galaxies. In the previous sections, we mea- 
sured the clump properties and carefully assessed the system- 
atic uncertainties associated with their extraction, in particu- 
lar those related to the treatment of the host galaxy's back- 
ground light. In this section, we now compare the properties 
of clumps in our NIC2 sample with those identified in other 
galaxy samples and at other wavelengths. We further explore 
constraints provided by our results on the origin of clumps 
and their role in the evolution of the host galaxies. 

5.1. Ensemble Properties of Clumps and their Origin 
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5.1.1. Comparison to Clumps in Other High-Redshift Galaxy 
Samples 

The properties derived in § |3] for the clumps in our 
SINS NIC2 targets agree well with the ranges for clumps 
in high-redshift galaxies derived from high resolution obser- 
vations at other wavelengths. The rest-frame optical sizes 
of ^ 1 kpc of our NIC2-identified clumps are comparable 
to those measured in the rest-frame UV and Ha line emis- 
sion of clumps identified in HST/ACS optical imaging and 
in AO-assisted SINFONI observations of z ^ 1-3 galax- 
ies (e.g., lE lmegreen & Elmegreen' 12005 at lElmegreen et al. I 
I2009at iGenzeletal. 2008. 2011). The fractional contribu- 
tions to the observed //leo band emission in our SINS NIC2 
objects, as well as to the observed /gi4 band emission in 
MD41 and Ha light in BX482, are typically a few percent 
per clump (though with a wide range of ^ 0.5% - 16%), 
and together make up ~ 10% -25% in each galaxy. These 
fractions are also very similar to the fractions in ACS ijjs 
band emission of res t-UV clumps in "clump- c lusters" and 
"chains," for which [ Elmegreen & Eknegreenl (12005 al) and 
lElmegreen et al.l (l20C)9al) derive tvpicallv ^ 2% per clump 
and a total of ~ 25% on average (ranging from < 1% up to 
^ 50%) for all clumps in a given galaxy. Clump-clusters 
and chains are the most prominently clumpy systems in op- 
tic al imaging of high redshift galaxies in the classification 
of lElmegreen. Elmegreen. & HirstI (I2004al) : rest-frame UV 
clumps in more regular galaxies at similar redshifts tend to 
have lower fractional contributions (Elm egreen et al. |[2009al) . 
In Ha emission, Genzel et al. (2011) determined compara- 
ble fractions in the range < 1% up to 25% for clumps in four 
z^2 disks. Thus, even if the clumps are a qualitatively promi- 
nent feature of the morphologies, the apparently smoother and 
more extended emission component generally dominates the 
total galaxy light. 

In terms of stellar m asses, Elmegr een & Elmegreenl ( l2005al) 
and lElmegreen et al. I ( 12009 a ) found that clumps in clump- 
clusters and chains appear to span a roughly constant range at 
z > 0.5 of MJ> ~ 10*^-5 X 10*^ Mq. The clump stellar masses 
in our SINS NIC 2 sample tend to l i e in t he top half of the mass 
distributions of Elmegre en et al. I (L2009a). This may reflect 
higher host stellar masses for most of our SINS NIC 2 targets 
compared to the gala xies studied bv|Elmegreen et al.l (see also 
the discussion bv iGenzel et al.ll2008l) . It mav also possiblv 
result from our clump selection in the rest-frame opti cal i n- 
stead of rest-frame UV as suggested by the analysis in § |4.1| of 
MD41, where NIC2-identified clumps are overall more mas- 
sive than ACS-identified clumps. Obviously, the important 
uncertainties involved in estimating the clump properties, and 
the differences in methodology — notably in the treatment 
of the background light, and in the specifics of the modeling 
used to estimate masses (see Elmegreen et al. 2009a, for de- 
tails) — mean that only order-of-magnitude comparisons can 
be made between the different clump samples. Nevertheless, 
these various studies indicate that individual giant clumps typ- 
ically represent each a few percent of the total stellar mass of 
the host galaxies. 

In the strong ly-lensed z ^ 1 -5 - 3 s tar-forming galaxies 
studied bv Jones et al.l (120101) and Hw inbank et al. (20103 
see also lA.M. Swinba nk et al. 201 1. in preparatiom). the sizes 
and dynamical masses of clumps identified in Ha or in sub- 
millimeter dust continuum and CO line emission tend to be 
lower (by factors of 2- 10) than for the clumps in our 
SINS NIC2 objects. Aided by linear magnification factors 



of 2-35, the observations for these lensed sources reach a 
source-plane spatial resolution down to w 100 pc in the best 
cases, undoubtedly helping to resolve smaller clumps. How- 
ever, the lensed galaxies are also on average roughly twice 
smaller and an order of magnitude less massive than the SINS 
NIC2 galaxies, and have lower inferred Vc/cto ratios (imply- 
ing dynamically hotter disks). The differences in clump prop- 
erties between the SINS NIC2 and the lensed samples may 
thus be more fundamentally related to the differences in host 
galaxy properties, as expected in the theoretical framework of 
gas-rich Toomre-unstable disks discussed next. 

5.1.2. Comparison to Theoretical Expectations for Fragmenting 
Gas-Rich, Turbulent Disks 

Theoretical arguments as well as numerical simulations of 
turbulent gas-rich Toomre-unstable disks indicate that they 
can fragment into large and massive star-forming clumps 
(e.g., iNoguchi 1999; Immeli et al. 2004a b; Boumaudetal. 
'2007; Genzel etal. 2008; Dekel et al. 2009; Ageitz et al 
| 2009; .Ceverino et al. 2010; Aumer et al...2010; .Genel et a" 
120101) . All our disks have high inferred local intrinsic gas 
velocity dispersions (^ 40-90 kms~'; iGenzel et al. I I2008t 
iCresci et al. 2009) and appear to be very gas-rich (gas-to- 
baryonic mass fractions of ~ 30% - 70%: lErb et al.l I2006t 
iForster Schreiber et al. I l2009t iTacconi et al. 112010b . Given 
these global properties, a possible interpretation of the clumps 
seen in our five disks is that they are the result of disk insta- 
bilitie s. Based on high S/ N AO-assisted SINFONI observa- 
tions, IGenzel et al.l (1201 II) argued that giant clumps seen in 
Ha and in rest-UV or optical broad-band emission of several 
large z ^ 2 star-forming disks (including BX 482) generally 
correspond to highly unstable regions in the galaxies. Indeed, 
IGenzel et al. I ( 1201 II) estimated from the AO data Toomre Q 
parameter values that are well below unity at the clump loca- 
tions. Moreover, Genzel et al. also inferred Q values in the in- 
terclump regions below unity, suggesting global perturbations 
across the host disks. These findings thus appear to be con- 
sistent with the hypothesis that the observed clumps formed 
in-situ from gravitational instabilities in gas-rich disks. 

In this framework, the characteristic size and mass of the 
forming clumps correspond to those of the fastest growing 
Jeans-unstable fragmentation mode not stabilized by rota- 
tion, the Toomre scale length and mass. These characteristic 
scales are related to global properties of the host galaxy and 
can be expressed as a function of disk radius and mass, cir- 
cular velocity Vd, and local intrinsic gas velocity dispersion 
(To (e.g., Genzel etal. 2008, 2011; Escala & Larson 2008] 
Elmegreen 2009bt iDekelet al. 11200 9). The ratio of va/ao 
is itself directly related to the disk scale length and scale 
height, and is inversely pro portional to th e gas mass fraction. 
Applying Equations (5) of IGenzel et al.l (1201 1) with Q ~ I 
and the assumption of constant rotation velocity, the global 
galaxy parameters for the five disks among our NIC2 sample 
(from Genzel et al. 2008; iCresci et al. 2009) imply a typical 
Toomre radius of ~ 1 kpc and mass of^4x 10^0. Compa- 
rable values were obtained bv I Genzel et al. I (12008b for a sam- 
ple of eight galaxies drawn from the SINS survey, with four 
objects in common with our NIC2 sample (MD41, BX389, 
BX610, and BX482). For Q < 1, these size and mass es- 
timates would increase (all other galaxy parameters being 
fixed). 

The Toomre length and mass inferred above are of the same 
order as the intrinsic sizes and stellar masses for the resolved 
clumps in the disks of our NIC2 sample. The total masses 
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could be significantly larger than the stellar masses if the 
clumps were themselves very gas-rich. Assuming that the 
typical gas-to-baryonic mass fraction of ^ 50% for the SINS 
N1C2 disks is representative for the clumps, the inferred to- 
tal mass for the resolved disk clumps would be in even closer 
agreement with the expectations. Direct molecular gas mass 
measurements for individual clumps in our targets are not 
yet available. However, estimates can be made for BX 482 
with AO-assisted SlNFONl data from Ha-based star forma- 
tion rates and using the "Kennicutt-Schmidt" relation. Fol- 
lowing Genzel et al. (201 1), to which we refer for details of 
the calculation and the choice of Kennicutt-Schmidt relation 
(slope and zero point) that might be most appropriate on the 
1 kpc-scales of clumps in high-redshift star-forming disks, the 
"background-subtracted" Ha measurements imply M^!,j, from 

< 6 X 1 0*^ Mp, for clump #2 to - 5 x 10'' Mq for clump #1 
(see also lGenzel et al. 1120 111 their clump "BX482-A"), with a 
median of ~ 10^ Mq (including limits). These estimates in- 
clude a 36% cor rection for helium and metals, account for the 
IChabried (120031) IMF adopted in this paper, and assume that 
the global Ay = 0.8 mag from the SED modeling of BX 482 is 
representative of the extinction towards the clumps. The gas 
mass for clump #1, the brightest amon g our N1C2 s ample, is 
comparable to those derived by Genzel et al. I (1201 Ih for other 
bright clumps in three additional z ^ 2 disks. For the N1C2- 
identified clumps in BX482, our results give Mg!,j. roughly 
twice the Mj' estimates, and total (stellar + gas) clump masses 
with median ^ 1.5 x 10^ Mq. The implied gas-to-baryonic 
mass fractions of ~ 60 - 70% for the clumps are comparabl e 
to the galaxy-integrated value dForster Schreiber et al. Il2009h . 

The theoretical expectations computed here as well as the 
empirical measurements should be regarded as rough esti- 
mates in view of the simplifying assumptions and the large un- 
certainties associated with the determination of clump proper- 
ties. In addition, differences between expectations and mea- 
surements can plausibly arise from different global galaxy 
properties at the time when the observed stellar clumps were 
formed and from evolution in clump properties (e.g., through 
mass loss). Notwithstanding these limitations, the agreement 
between the predictions and the observed properties is re- 
markable, and provides a framework to interpret the promi- 
nent clumpy feat ures in distant gala xies. Similar conclusions 
were reached by I Jones et aP (l201Cft for the origin of clumps 
in their lensed galaxy sample, based on the measured clump 
sizes and the expectations computed with the size, mass, and 
dynamical prope rties of the host d isks. As pointed out in 
previous studies dfescala & LarsonI (20081: lElmegreenI l2009bl: 
iGenzel et al. 11201 IL and references therein), the high gas mass 
fractions and low Vd/ao ratios observed in z > 1 disks nat- 
urally lead to the formation of larger self-gravitating, star- 
forming complexes than in the less gas-rich, dynamically 
cold, and geometrically thin z ^ disks. In addition, the 
presence of a massive stabilizing stellar disk and/or bulge 
in local mature disk galaxies also leads to smaller lower- 
mass star-forming comp lexes, or even to the suppression 
of their formation (^e.g.. iBournaud et al. I l2007t iDekel et aP 
l2009l:lCeverino et al.ll20m 

5.2. Constraints on Radial Age Variations and Clump 
Evolution 

In the scenario outlined in the previous section, the mas- 
sive star-forming clumps are expected to migrate towards the 
gravitational center as a result of dynamical friction against 



the host disk and of clump-clump interactions. Ultimately, 
the clumps may coal esce into a young bulge within ^ 10 
dynamical timescales ( Noguchilll99 9'; 'Immeli et al. "2004af3; 
Boumaud et al. 200 71: iGenzel et al. . 2008: Dekel et al. 20091 
Agert z et al. I 120091: ICeverino et al. 1 1201 0) unless they are 
rapidly disrupted, for instance by stellar feedback or tidal 
torques. In that respect, recent theoretical work has high- 
lighted the possibly dramatic impact on clum p evolution 
of vigorous star formati on feedback (M urray e t al. 1 120 10[ 
iGenel et aDl2010l but see Kru mholz & Deke]ii201C f for a con- 
trasting view). High resolution and high S/N SINFONI ob- 
servations have provided the first empirical evidence of gas 
outflo ws originating frorn massive luminous clumps in z 2 
disks (IGenzel et al. 11201 lb . In some of the most actively star- 
forming clumps, the outflow could be sufficiently strong to 
disperse a large fraction of the initial gas before they reach 
the center of the galaxies. If the clumps do survive long 
enough and spiral inward, they are then expected to exhibit 
a significant age spread, with olde r clumps typically closer 
to the center of the gal axy (e.g., lElmegreen et al. I l2009at 
iKrumholz & Deke]||2010l) . 

Our results for BX482 seem broadly consistent with this 
picture. From measuremen ts of the age-sensitive Ha equiv- 
alent width, we found in § 14.21 indications of age variations 
among the NIC2- and Ha-identified clumps. Figure |8] shows 
the distribution of W^'^^'(Ha) and of the corresponding stellar 
ages assuming CSF as a function of deprojected galactocen- 
tric distance d. Although the data do not reveal a significant 
correlation between W'^^'(Ha) or age and d (a Spearman rank 
correlation analysis yields a coefficient of about 0.65 but with 
a significance of « 1 .5 a), there is a clear distinction dri ven by 
the NIC2-identified clump #4 already highlighted in S I4.2H 
This clump, with the lowest Ha equivalent width and oldest 
age, is also nearest the center of BX482 at d = 2.9 kpc. The 
six other clumps, with higher equivalent widths and younger 
ages, are located at significantly larger d « 5-7.5 kpc. This 
distinction is qualitatively preserved for alternative choices 
of SFH as long as the clumps follow similar SFHs (see Ap- 
pendixlAli. 

For MD41, the stellar ages cannot be constrained directly 
but the deprojected radial distribution of isi4-H[(,o colors of 
the clumps can provide indications of maturity in terms of 
M/L ratio. The colors are plotted as a function of galacto- 
centric distance d in Figure|9] along with the derived M,t/L^''^' 
ratios. These plots suggest a weak trend (with Spearman rank 
correlation coefficient of « 0.85 and significance of 2 a) of 
redder colors and higher M/L ratios at smaller radii; it is 
mostly do minated by clump NIC2-#1 and ACS-#3, which, as 
noted in § 14.1.21 exhibit the most distinct colors compared 
to the ensemble of clumps and interclump regions B If the 
dust extinction and the SFHs are similar among the clumps, 
the radial distribution of colors and M/L could be indicative 
of older stellar ages for clumps at smaller galactocentric dis- 
tances. 

In their studies of clumpy z ^ 1 -4 galaxies in the H-UDF, 
lElmegreen & ElmegreenI (l2005al from ACS B435, Veoe, '775, 

The distribution of W^HHa) with d over the interclump pixels in 
BX482 shows a large scatter with hardly any trend (p w 0.18 at the 2(T level), 
reflecting the over all ve ry similar appearance of the Ho and //igo band emis- 
sion discussed in § 14.21 

The distribution of igi4 —H[^q colors with d over the interclump pixels 
in MD41 is almost flat, with increasing scatter at larger radii (and p 0.15 
at the 2 a level), as indicated by the fairly uniform colo rs ac ross MD 4 1 with 
the exception of some of the clumps as highlighted in § 14.11 
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and zs5o band imaging) and lElmegreen et al.l (l2009ai using 
the same ACS bands, complemented with NICMOS/NIC3 
7iio and i/ieo band imaging) found that ACS 1775 band- 
identified rest-frame UV clumps within a given galaxy mostly 
span a narrow range of colors and that the interclump regions 
tend to be redder, which they interpreted in terms of older ages 
for the underlying stellar population across the host galaxy. 
In ~ 30% -50% of the H-U DF chains and clu mp-clusters, 
and in all their "spiral" types, lElmegreen et al. I (^009a) also 
identified clumps that are prominent in the near-lR bands, 
exhibit redder colors, and tend to lie closer to the galaxies' 
centers. These "bulge-like clumps" (or bulges for spirals) 
have higher inferred stellar ages and masses than the other 
clumps present within the same galaxy; they do not appear 
to differentiate significantly in terms of star formation decay 
timescales or extinction. Moreover, these differences in stel- 
lar ages and masses are more important for spirals, less signif- 
icant for clump-clusters and chains. The derived stellar ages 
of clumps span a wide range, center ed aroun d ^100 Myr. As 
argued bv iElmegreen et aO (l2009ai see also Elmegreen et al. I 
l2009cl) . the inferred clump ages and the presence of a red- 
der and more massive clump generally closer to the center in 
clumpy and spiral types are consistent with inward clump mi- 
gration, and suggest that at least a fraction of the clumps may 
survive long enough to reach the central few kpc. 

The trends of redder colors for the clumps in MD41 and 
of older ages for those in BX482 at smaller galactocentric 
radii are consistent with the general findings of Elmegreen, 
Elmegreen, and coworkers. Our higher resolution NIC2 imag- 
ing, compared to N1C3, allowed us however to detect clumps 
directly from longer wavelengths on ^ 1 kpc scales. The 
N1C2 data furthermore provide a better match to the ACS 
PSF enabling us, in the case of MD41, to identify and mea- 
sure the properties of clumps consistently in both the opti- 
cal and near-lR bands. While most ACS -identified clumps 
in MD41 tend to be bl uer than the int e rclump regions, in 
Une with the results of .Elmegreen et al.l (l2009al) . the N1C2- 
identified clumps have in contrast comparable or redder col- 
ors. The small but systematic differences in /8i4-i^i60 col- 
ors, and in derived M/L ratios and stellar masses, between the 
NIC2- and ACS-identified clumps would suggest that clumps 
identified in the rest-optical and rest-UV are complementary 
and probe together a wider range in properties. We can 
draw a similar conclusion for the clumps identified in rest- 
frame optical and Ha in BX482, with the N1C2 and SIN- 
FONl PSFs also well matched to each other, although there 
is more overlap here. The clump ages in BX482 for CSF 
models are compatible \y ith the range for clumps at z ~ 2 by 
lElmepreen et al. l (l2009al) : clump NIC2-#4 closest to the cen- 
ter lies at t he high age end while all others are in the younger 
half of the lElmegreen et al.l distribution. 

The dynamical timescale t^yn = r/v at the effective radius 
w 6 kpc of BX482 is w 25 Myr (correspondingly, the or- 
bital timescale is rorb('" = ^eff) ~ 160 Myr). If the clumps have 
been forming stars at a roughly constant rate, the stellar age 
of the oldest clump NlC2-#4 thus implies it has survived for 
~ 10 to ^20 dynamical timescales, depending on whether 
the raw or background-subtracted Ha equivalent width esti- 
mates is adopted. It thus appears plausible that this clump has 
formed at larger radii and has migrated to its present location 
closer to the center at li = 2.9 kpc. The other younger clumps 
have ages ~ (1 -4) x f^yn, so that they may still be located 
fairly close to the radius at which they formed. For the largest 
and brightest clump NIC2-#l. ,Genzel et al..(.2011.) detected a 



broad (FWHM w 400 kms"') and somewhat blueshifted (by 
several 10 kms~') Ha Une component, which is interpreted as 
a sign of outflowing gas. Combining estimates of the outflow 
kinematics and rate, and of the clump extent, gas mass, and 
gas-phase oxygen abundance, iGenzel et al. I 02O1 II) derived 
gas dispersion, outflow expansion, and chemical enrichment 
timescales in the range ~ 0.3- 1 Gyr. These estimates would 
imply longer clump lifetimes than our W^'''*'(Ha)-based stellar 
age for CSF of « 55 Myr, and may suggest that this clump has 
survived for more than a couple dynamical timescales. 

5.3. Constraints on a Stellar Bulge/Inner Disk Component 

Clump migration may provide an important mechanism for 
growing a stellar bulg e and/or inner disk component in hig h- 
redshift galaxies (e.g. 'Noguc hill 19991: llm meli et al. "2004a'lbl: 
Boumaudetal. 2007; Caroll o et al. 1 120 07; Elmegreen et alT] 
2009a; Ce verino et al.l 120 lOi) . along with other internal 
mechanisms such as stellar bar-induced streaming (e.g., 
iBourna ud & Combes 200 3) or visco us drag in a gaseous disk 
(e.g.. iLin & Pringl&.1987Jsilkll2001i) . In z ~ disks, these 
secular pro cesses occur on long timescale s of at least several 
Gyrs (e.g., 'Kormend y & Kennicut3l2004 '). At z - 2, how- 
ever, disk galaxies are significantly more gas-rich, and their 
gas phase appears to be much more turbulent. Under these 
conditions, secular processes should proceed faster by an or- 
der of m agnitude or more (e.g. Noguchi 1999; Immeli et aQ 
'2004a','b'; 'Boumaud etJQ l2007c iGenzel etal.i ,2008l ,201 It 
Ceverino et al. 2010). If clumps are dispersed and lose mass 
before reaching the very center of the galaxy, the build-up 
of a bulge component through clump migration would be 
slower but this mec hanism could still contribute to inward 
mass transfer (e.g., iBoumaud et al. 1 120071: lElmegreen et al. I 
120093) . 

From duty cycle arguments, lElmegreen et al. I (l2009al) in- 
ferred that the clumpy phase in high-redshift disks may last 
for several Gyrs, or 5-10 clump formation epochs. Some 
recent high-resolution cosmological simulations suggest that 
rapid accretion from the halo along narrow, cold streams re- 
plenishing the gas reservoir of high-redshift disks could help 
to maintain the clumpy phase over such extended timescales 
(iDekel et al.ll2009tlCeverino et al.l lMo^. One may therefore 
expect the presence of a more important stellar bulge or in- 
ner disk component in more mature galaxies, as inward mi- 
gration of previous clump generations and secular processes 
would have had more time to transport material in the central 
regions. 

In this context, an immediate question is whether our data 
provide evidence for a more important stellar bulge or inner 
disk component in our NIC2 sample disks with older stel- 
lar populations. At first glance, such a trend is not appar- 
ent from the observed spatial distribution of the rest-frame 
optical continuum emission seen in our NIC2 //i6o-band im- 
ages. The best-fit Sersic index n does not show any significant 
correlation with best-fit stellar age from the rest-UV/optical 
SED modeling, as discussed in iPaper 3 . The implied stellar 
mass fraction of the clump closest to the center in each of the 
disks, if interpreted as a young bulge, also shows no trend 
with galaxy stellar age. More generally, none of the disks 
in our sample shows obvious signs of a prominent centrally- 
concentrated stellar component apart fro m BX 66 3, for which 
the two-component Sersic model fits in iPaper II suggest that 
the central peak could be associated with a fairly compact 
« = 4 bulge-like component contributing w 35% of the total 
//i60-band luminosity. However, there could be an older stel- 
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lar population in the inner parts of the galaxies eluding detec- 
tion even at rest-frame « 5000 A if it has a significantly higher 
M/L ratio than the stellar population dominating the observed 
rest-frame UV to optical SED. 

To set constraints on the mass contribution of such a hid- 
den stellar component, we postulated that this component 
has an effective radius = I kpc and Sersic index n = 3 
akin to the compact quiescent z ~ 2 /T-selected ga laxies of 
lyan Dokkum et al. ( 2008, see also lKriek et al.ll2009h . and as- 
sumed a fixed axis ratio b/a= 1. In what follows, we refer to 
this component as "bulge," motivated by our choice of struc- 
tural parameters, but we note that other interpretations are not 
excluded. We estimated how much extra light from this com- 
ponent is allowed by the 1 cr uncertainties of the observed ra- 
dial light profiles at r < 2 kpc for each of the five disks among 
our N1C2 sample, as illustrated in Figure [TO] To determine 
the corresponding stellar mass contributions, we considered 
in turn different M^/L^'^'*' ratios. Arguably, a realistic choice 
is the t ypical M/L found for z ~ 2 compact quiescent galax- 
ies (see lKriek et aL]l2009l) . For the Ivan Dokkum et ail (120081) 
sample, the medi an Mj,/Lg'^' [ MQL~ '(q] sa 0.6 based on the 
SED mod eling of lKriek et al.l ( l2008h . after adjustment to the 
iChabried (120031) IMF adopted in this workEJ This value is 
uncorrected for dust extinction, but the median Ay for the qui- 
escent sample is similar to that of our five SINS NIC2 disks 
(0.7 versus 0.8 mag, respectively). This "passive" M/L is, 
however, comparable to or lower than the ratio for our three 
disks with oldest stellar ages. If there was an important bulge 
characterized by such a M/L ratio, it would be clearly seen 
in the data. We therefore explored the more extreme scenar- 
ios in which all stars of the hidden bulge are as old as the 
age of the universe at the redshift of each disk, and whose 
light is either unobscured or attenuated by the same amount 
as the stars that dominate the rest-UV/optical SED. The M/L 
for this maximally-old single stellar population (SSP) is cal- 
culated from the same Bruzual & Chariot (2003) models for 
solar metallicity as used in Paper I for the SED fitting. 

The results are reported in Table |5] The light contribution 
7^''"'s'^(//if,o) corresponds to the fraction of the total observed 
i/i6o-band flux density for each galaxy, and the mass fraction 
is relative to the stellar mass obtained from the SED model- 
ing, J"'"''8'=(M,,) = M^'*^7Mf D. The uncertainties given in 
Table |5] reflect the (small) impact of varying the bulge Ser- 
sic index between « = 2 and 4, and of increasing the up 
to 3 kpc. Unsurprisingly, the light fractions are very small, 
in the range l%-3.5%. In terms of stellar mass, the highest 
fractions hidden in an older bulge are found for the disks with 
youngest ages (hence lowest Mi,/U^^^ ratios) from the best- 
fit SED model. Nevertheless, the mass fractions are small to 
moderate only, ranging from a few percent to w 30%, with the 
exception of MD 41 (with estimates of up to a factor of ^ 2.3, 

but see below). For all objects, the sur n of M^^°+M^"'^^ 
does not exceed the total dynamical mass (iGenzelet al.ll2008l; 
ICresci et al.|[2009h . even whe n accounting for the substan- 
tial amounts of gas infeiTed (iForster Schreiber et al. 1 120091 
iTacconi et al.ll2010l) . A further plausibiUty check for our as- 
sumptions can be made for MD 41 by comparing the resulting 
M/L ratio for the various cases with the /814-//16O color-based 
estimates in the inner r < 1 kpc. Only the maximally-old ob- 

^'^ The difference between M^/U^^ relevant for our data and M*/LJ^^' 
given bv lKriek et all is < 10% for the typical ages, Ay, and SFHs of 

the quiescent galaxies. 



scured bulge scenario leads to a Mi,/Ug^^ value that is largely 
different from the ratio corresponding to the observed col- 
ors, suggesting that this case (with more than twice the stellar 
mass being hidden) may be regarded as unlikely. More gen- 
erally, the maximally-old SSP assumption is an extreme one 
and the results should be considered as upper limits. 

We conclude from the exercise above that the disks in 
our NIC2 sample are unlikely to contain an important frac- 
tion of their stellar mass in an old bulge or inner disk un- 
seen in our //leo-band data, with the possible exception of 
MD41. BX663 appears to be the only object in which a sig- 
nificant bulge-like component is directly visible, if the AGN 
contribution to the rest-frame k, 5000 A emission can be ne- 
glected. The visible and hidden mass estimates add up to 
~ 1 X 10''-3 X 10'° M0 in the centi'al r < 1 kpc regions of 
our SINS NIC2 disks; thus, it appears that none of these disks 
hosts a compact massive core. No trend between stellar pop- 
ulation age and the importance of a central stellar component 
emerges when a ccounting for a n old hidden population. 

Interestingly, Genzel et al. I (12008 ) found in contrast evi- 
dence for a connection between central dynamical mass con- 
centration and maturity of the stellar population, with the ratio 
of dynamical mass enclosed within r « 3 and 10 kpc increas- 
ing with older stellar age and higher inferred metallicity of the 
galaxies. For the four objects in common with our NIC2 sam- 
ple, Genzel et al. derived Mdyn('" < 3kpc)/Mdyn(r < lOkpc) < 
0.15 for MD41, 0.21 for BX482, and 0.39 for both BX389 
and BX610. Combining the observed //leo-band light pro- 
files from our NIC2 data together with the constraints on the 
mass of a hidden evolved bulge, we find that the implied stel- 
lar mass concentration calculated in the same manner agrees 
reasonably well with the value inferred from the dynamics for 
BX482 and BX610. On the other hand, there are significant 
differences between the central stellar and dynamical mass 
concentrations for MD 41 and BX 389. 

Thus, despite our findings in § |4] consistent with inward 
clump migration — a possible mechanism for growing stellar 
bulges or inner disks over time — our SINS NIC2 disks do 
not appear to show a trend of higher central stellar mass con- 
centration with older galaxy age even when accounting for a 
possible hidden component. Admittedly, our sample may be 
too small, and may sample too sparsely different evolutionary 
stages between the objects with younger 50-320 Myr) 
and older 2.5-2.8 Gyr) stellar ages, to discern such a 
trend. Moreover, various factors could obviously complicate 
the analysis presented in this section, including the unknown 
spatial variations in extinction and in the distribution of the 
gas, which contributes a large fraction of the total mass in 
these galaxies. Detailed maps constraining the extinction and 
gas mass distribution, in addition to the M/L ratio and the 
kinematics, will be important to further test the scenario of 
bulge/inner disk growth through clump migration and secular 
processes. 

5.4. Other Interpretations of the Clumpy Features 

In this paper, we assumed that the compact regions of 
higher apparent rest-frame optical surface brightness across 
our SINS NIC2 targets represent genuine substructure in the 
distribution of the stellar populations. As we discussed earlier 
in this section, the properties derived for these clumps as well 
as the global properties of the host galaxies are consistent with 
the scenario in which giant massive clumps result from gravi- 
tational instabilities in gas-rich, turbulent disks. This scenario 



14 



Forster Schreiber et al. 



provides a framework to interpret clumpy morphologies and 
to investigate dynamical processes that may play an important 
role in the build-up of bulges and inner disks in high-redshift 
galaxies. 

Other interpretations are possible. Bright clumps could rep- 
resent merger components. BX 528 is a case where at least 
two of the clumps identified most likely correspond to the 
progenitor galaxies of a major merging pair, and the southern 
clump # -4 in BX 389 is plausibly a small companion galaxy 
(iPaper HI) . Clumps within disk galaxies could also be of exter- 
nal origin, and correspond to newly accreted low-mass sys- 
tems. The distinction between clumps formed in-situ and 
those accreted from the halo is not easy. Constraints could 
be obtained from an accurate analysis of the stellar proper- 
ties of clumps and host galaxies including age and metal- 
licity, or from a detailed baryonic and dynamical mass bud- 
get of clumps to test whether they contain dark matter Ar- 
guably, it may seem however rather unlikely that as many as 
5-10 clumps per galaxy all originate from rapid accretion of 
satellites or clumps brought in from the halo along filaments 
or dense, cold narrow streams (e.g ., Elmegreen et al. 2008; 
iDekel et al.ll2009tlAgertz et al.ll2009l:iCeverino et al.il2010i) . 

Another possibility is that the kpc-sized clumps simply cor- 
respond to locations across galaxies with lower line-of-sight 
dust obscuration. Here again, the distinction between real 
substructure in the stellar populations or star-forming regions 
and regions of lower obscuration would necessitate detailed 
spectral or multi-wavelength constraints on resolved scales of 
1 kpc or smaller to pin-down accurately the variations in ex- 
tinction across the galaxies. The role of extinction in clumpy 
morphologies is certainly a more important concern in stud- 
ies based on rest-frame UV imaging compared to rest-frame 
optical imaging, but it could still play a significant role at 
rest-frame « 5000 A, as is well-known from nearby dusty 
star-forming galaxies. Small fractions of light of a few per- 
cent, as inferred for most clumps in our SINS N1C2 sam- 
ple as well as on average for clumps in optical imaging of 
high-redshift galaxies, do not rule out this interpretation, al- 
though it appears less plausible in the case of very bright 
clumps with fractions in excess of ^ 10%. Extinction effects 
may also naturally explain the findings by lElmegreen et aP 
( l2009a) that rest-UV identified clumps tend to have bluer 
colors than the interclump regions. It may be more diffi- 
cult to reconcile with the trend of redder colors we found for 
rest-optically identified clumps in MD41, or the "bulge-like" 
cl umps often at more centra l locations in galaxies discussed 
by lElmegreen et aD (l2009a l). If large localized spatial varia- 
tions in extinction are responsible for the clumpy features in 
our galaxies, trends of redder colors and old er ag es at smaller 
radii, as inferred for MD 41 and BX 482 in § |42] could reflect 
age gradients in the underlying stellar population, dominated 
by younger stars at larger radii and older ones in a bulge or 
inner disk component. 

Possibly, clumpy features in high-redshift galaxies have 
multiple causes. In future studies, high spatial resolution 
multiwavelength mapping and sensitive spectroscopy will be 
essential to determine robustly the nature of clumps. More 
realistically in the near future, a combination of optical and 
near-lR imaging (allowing spatially-resolved SED model- 
ing), AO-assisted integral field spectroscopy (mapping the 
kinematics and star formation, and providing additional con- 
straints on the age and SFH), and millimeter interferometry 
(mapping directly the molecular gas distribution and kine- 



matics) will allow tighter constraints on the detailed physical, 
stellar, and dynamical properties o f clumps. The first steps 
taken in this direction (e.g., lElme green & ElmegreenI 
2005a; Elmegreen etaP l2009d iGenzel et al. 20liF 
Tacconi et al. 2010; Jones et al. 2010; Swinbank et aD 
2010b ; A.M. Swinbank et al. 201 1, in preparation) show 
that, if challenging, such an approach will provide key 
information necessary to better understand the nature of the 
observed kpc-scale substructure and its role in the evolution 
of distant galaxies. 

6. SUMMARY 

We have studied the characteristics of the clumpy features 
based on HST/NIC2 //leo band imaging for a unique sample 
of six z ^ 2 star-forming galaxies with VLT/SINFONI integral 
field spectroscopy. The N1C2 imaging provides good sam- 
pling of the HST PSF at 1.6 /im and allows us to detect reli- 
ably substructure in the rest-frame optical emission on scales 
as small as 1.2 kpc at the z = 2.1 -2.5 of our targets. We 
have identified "clumps" through a systematic detection pro- 
cedure that takes into account both host galaxy background 
subtraction and PSF deconvolution for estimating the clump 
luminosities, stellar masses, and sizes. We assessed for the 
first time in detail the systematic uncertainties associated with 
extracting clump properties in distant galaxies, especially in 
terms of the treatment of the underlying host galaxy light. We 
interpreted the properties of the clumps in the framework of 
gas-rich, turbulent disks in which giant massive clumps can 
form through violent gravitational instabilities. Our study 
complements previous work focussing on clumps identified 
in the rest-frame UV or in Ha line emission, and shows the 
importance of including constraints at longer wavelengths and 
on the same ^ 1 kpc scales for a more complete picture. 

Our main results can be summarized as follows. 

1 . We identified between 2 and 7 clumps per galaxy in their 
rest-frame optical emission. Using the SlNFONl data, we 
determined that these clumps are unlikely to be caused pre- 
dominantly by strong nebular emission lines within the //leo 
bandpass and thus must reflect compact sources of stellar con- 
tinuum light. At least two of the clumps in the kinematically- 
identified binary merger BX 528 are likely associated with the 
merging galaxies. For the large kinematically-confirmed disk 
BX389, the clump most offset from the main body of the 
galaxy can be securely identified as a small nearby compan- 
ion at a projected distance of « 5 kpc. Continuum light from a 
stellar bulge component could plausibly dominate the central 
prominent peak in BX663. In terms of size, luminosity, and 
stellar mass, these specific clumpy features do not differen- 
tiate significantly from the ensemble of all clumps identified 
in our N1C2 targets, emphasizing the importance of taking 
into account the nature of the host system from the dynamics, 
morphologies, and spectral properties in the interpretation of 
"clumps." All other clumps may have plausibly formed in- 
situ from fragmentation of the turbulent, gas-rich disks of our 
targets. 

2. The stellar clumps in the //i(,o band images of our N1C2 
targets generally represent surface brightness enhancements 
of ^ 25% - 35% relative to the surrounding background from 
the host galaxy. This contrast is substantially lower than for 
clumps identified in optical ACS imaging of z ~ 1—4 promi- 
nently clumpy galaxies (clump-clusters and chains), which 
generally appea r to be brighter than the b ackground by factors 
of ^ 2-4 (e.g.. Elmegreen et al.ll2009"al) . The contribution of 
individual clumps to the total galaxy i/igo band emission is 
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of a few percent typically, although with wide range from the 
faintest (w 0.5%) to the brightest (w 16%). The total contri- 
bution of all clumps in a given galaxy is between 10% and 
25%. The typical intrinsic rest-frame optical FWHM size of 
our clumps is ^ 1 kpc, similar to sizes of clumps identified 
and measured in the rest-UV conti nuum and Ha line emis- 
sion (e . g., [E lmegreen & Elmegreenl l2005at lElmegreen et al. I 
I2009at iGenzel et al., ,2008. ,201 ID . The inferred stellar 
masses have a median 10^ Mq, and a range that lies 
within the more massive half of the distribution inferred for 
rest-UV identified clumps (Elmegreen & Elmegreen 2005^ 
lElmegreen et al.ll2069al) . possibly reflecting the higher stellar 
masses of our galaxies and the clump identification at longer 
wavelengths. 

3. By combining the N1C2 data with H ST ACS ku band 
imaging available for MD41, and existing AO-assisted SIN- 
FONl Ha data for BX482, we infeiTed modest color, M/L, 
and stellar age variations among most clumps within each 
galaxy, and between clumps and interclump regions. In these 
two objects, the sets of clumps identified at different wave- 
lengths do not fully overlap. In MD41, NlC2-identified 
clumps tend to have redder /814-//16O colors than ACS- 
identified clumps and than the interclump regions. The col- 
ors for the reddest and bluest clumps appear to require dif- 
ferences in stellar age in addition to extinction effects. In 
BX482, we derived more robust constraints on the relative 
stellar ages from the Ha equivalent width, taken as the Ha 
line to rest-frame « 5000 A continuum ratio. Assuming sim- 
ilar star formation histories, the NlC2-identified clumps tend 
to be older than the Ha-identified clumps but, except for one, 
have ages younger than the interclump regions. These find- 
ings suggest that clumps identified at different wavelengths 
are complementary and probe together a wider range in clump 
properties. Most importantly, clumps in MD41 and BX482 
that are closer to the center appear to show more maturity in 
terms of stellar populations. 

4. The quantitative constraints on clump properties for the 
kinematically-confirmed disks in our SINS N1C2 sample are 
consistent with the scenario in which giant massive clumps 
form from disk instabilities. Their sizes and masses are of 
the order of the Toomre scale length and mass expected for 
disks with the same global properties as measured for the 
galaxies (large gas velocity dispersion, high gas mass frac- 
tions, and disk sizes). The trends of older and/or more ob- 
scured stellar populations of clumps at smaller galactocen- 
tric radii for BX482 and MD41 are consistent with inward 
migration of the massive clumps through dynamical friction 
and clump-clump interactions. These findings suggest that 
a fraction of the clumps can survive to destruction by stel- 
lar feedback or other disruptive mechanisms for at least a 
few dynamical timescales. Our results add to other evi- 
dence fro m previous studies o f high-redshift clumpy g alax- 
ies (e.g., IGenzel et al.l 120081: lElmegreen etaTl l2009al) and 
from num erical simulations of gas-rich uns table disks (e.g. , 
Noguchi|[T9 99; Immeli et al. 2004a b; Bour naud et aL]l2007l: 
Agertz et al. ii2009t iCeverino et al.i 12010 ) that clump migra- 
tion could be an important mechanism by which material is 
transported in the inner regions of young disk galaxies, con- 
tributing to the build-up of a bulge and inner disk component. 



5. We estimated for the disks in our sample ^ 1 x 10''- 
3 X 10'" Mq in a stellar bulge or inner disk component at 
r < 1 - 3 kpc, with no obvious trend emerging among these 
five sources between central stellar mass fraction and stellar 
age of the galaxy. Our galaxies appear to have roughly at least 
an order of magnitude lower stellar central densities than the 
compact quiescent objects and even some of the star-forming 
galaxies with significant putative compa ct cores among the 
massiv e 2 < z < 2.6 NIC2 s amples of Ivan Dokkum et alH 
(120081) and lKrieketaT] (120091) . This could suggest different 
path s for the forma tion of the inner ste llar components (see, 
e.g.. iTacconi et al.l [2008; Fra nx et all I2OO8, and references 
therein), with highly dissipative mechanisms (such as major 
mergers) being more important in the latter samples and sec- 
ular processes dominating for our SINS N1C2 sample. 

The sample presented in this work is small and spans a 
modest range in galaxy parameters, but our analysis high- 
lights some of the key issues that will be important to ad- 
dress in future studies. The identification of clumps at differ- 
ent wavelengths can potentially probe a wider range in clump 
properties and evolutionary stage. Larger samples probing a 
wider range in host galaxy properties and including (major) 
mergers as well will clearly allow to better explore trends of 
clump properties with galaxy parameters, evolutionary stage, 
and location within the host galaxy, and to assess possible 
differences between clumps in disks and in merging systems. 
The combination of M/L-sensitive and age-sensitive proper- 
ties will be crucial to make progress. One of the currently 
most important limitations in the interpretation of the clump 
properties is their unknown star formation history. Tighter 
constraints on clump masses, star formation and outflow rates, 
and internal structure and dynamics will provide better in- 
sights into the issue of clump lifetimes and their role in the 
internal dynamical evolution of young galaxies. 
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APPENDIX 

A. VARIATIONS OF CLUMP STELLAR AGES WITH DIFFERENT STAR FORMATION HISTORIES FOR BX 482 

In §§ l4.2l and l5.2l we used the Ha equivalent width, with the continuum flux density approximated from the i/ieo-band photom- 
etry of our targets, to derive stellar ages for the the clumps and interclump regions in BX482. This indicator is nearly insensitive 
to interstellar dust reddening if differential extinction is negligible, i.e., if the Ha line emission from the H II regions is attenuated 
by the same amount as the continuum light from the stellar population. On the other hand, the W^'''*'(Ha) depends sensitively on 
the star formation history (SFH). The analysis throughout the paper relies on the simple assumption of constant star formation 
(CSF) with age. Here we explore the effects of adopting other SF Hs on the interpretation of the W^^'C Ha). 

For this purpose, we u sed the same i Bruzual & Chariot! (120031) models of solar metallicity, with the I Chabrierl (12003 1) IMF and 
the lCalzetti et al.l (120001) reddening law as for the SED modehng of the integrated rest-UV/optical SED of our SINS NIC2 sample 
galaxies. We investigated three different generic forms for the SFH: (1) exponentially declining models with SFR(f) oc exp(-f/r) 
("t" models), (2) exponentially increasing models with SFR(f) oc exp(f/r) ("inverted r", or "ir" models), and (3) a hybrid case 
with SFR(f) oc (t /t) exp(-f /r) ("delayed r", or "dr" models), where SFR is the star formation rate, t is the time elapsed since 
the onset of star formation (taken as the age of the stellar population), and r is the characteristic timescale of the star formation. 
A CSF corresponds, for the exponentially-declining parametrization, to t — > oo. We computed suites of models for a range of r 
from 10 Myr to 1 Gyr For all SFHs, the W^'^^'(Ha) decreases monotonically with age; differences are in the shape of the curves. 

We derived the stellar ages and uncertainties of the clumps in BX 482 from the W"^^'(Ha) and their 1 a errors (projected 
onto the model curves). We always assumed a fixed maximum age possible corresponding to the age of the Universe at the 
redshift of BX482. Ages older than the Universe are formally obtained only for NIC2 clump #4, with lowest W^'^^'(Ha), and 
inverted r models with timescales > 100 Myr. This suggests that either such SFHs are inappropriate or there is significant 
differential attenuation between the Ha line and stellar continuum emission; the data currently available do not allow us to 
distinguish between these two possibilities. Figure [TT] illustrates the differences in ages obtained for each SFH family and 
three representative timescales, t = 30, 100, and 300 Myr, plotted as a function of the ages for CSF models. Because of the 
monotonic variations with t irrespective of the SFH, the relative age differences between clumps are qualitatively preserved for 
any given SFH. The range of ages tends to become narrower for shorter star formation timescales. However, the uncertainties 
of the measurements are typically sufficiently small that the age variations among the clumps remain significant for all r values 
considered. This result holds independently of the treatement of the background light from the host galaxy; the measurements 
based on the "background-subtracted" photometry lead to larger age differences than those based on the "raw" photometry, i.e., 
not accounting for the contribution of the underlying host stellar population to the clump light. 

From this exploration of the impact of SFHs on the interpretation of the W^'^^'(Ha), we conclude the following. While the abso- 
lute ages of the clumps in BX482 and ranges thereof depend on the SFH, the relative age variations are qualitatively robust and 
quantitatively significant in view of the measurements uncertainties. This conclusion is valid as long as the SFHs of individual 
clumps are similar in their general time dependence and in their timescale. Verifying this assumption would be important but 
obviously, detailed constraints on clump SFHs would require much more extensive multi wavelength data sets at high angular res- 
olution. Future studies with HSTAVFC3 and later with JWST and the ELTs should be able to provide such constraints and further 
insights in the detailed star formation history across individual galaxies, thereby also allowing more accurate reconstruction of 
the galaxy-wide global SFH, one of the major uncertainty in SED modeling of distant objects. 
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TABLE 1 

Galaxies observed 
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Kinematic Type'^ 


Major merger 


Disk 


Disk 


Disk 


Disk 
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0.56 ±0.01 


0.48 ±0.01 



" Spectral signatures of a Type 2 AGN are detected in the optical and near-I R spectnmi of BX 663 a nd the large-scale Ha kinematics 
suggest the host galaxy is a low-inclination rotating disk ( E rb et al. 20061 : iForster Schreiber et al. I t2006, 2009). The 5p!/zer/MIPS 
24 /im flux in BX663 is consistent with obscured AGN activity. A 24 /xm excess is also detected for BX610 but since there is no 
other signature of AGN activity in this galaxy, it could be due to enhanced mid-IR PAH emission around rest-frame 8 /im instead 
(see §|I2}. 

Systemic vacuum redshift derived from the integrated Ha line emission ( IForster Schreiber et al. Il2009h . 

The //i6o band magnitudes (AB system) measured from the NIC2 data in a circular aperture with diameter of 3". 

Stellar age, visual extinction (Av), stellar mass (M*), stellar mass to rest-frame SDSS g-band luminosity ratio (M*/Lj", uncor- 
rected for extinction), and absolute and specific star formation rate (SFR and sSFR) of the best-fit evolutionary syn thesis model 
to the observed optical to near-IR photometry of the galaxies. The SEDs were fitted using iBruzual & Charlo3l200 3' models with 
solar metallicity, a Chabrier 2003 IMF, the Calzetti et al. 2000 reddening law, and a constant SFR. The formal (random) fitting 
uncertainties listed are derived from the 68% confidence intervals based on 200 Monte Carlo simulations; systematic uncertainties 
(from SED modeling assumptions) are estimated to be typically a factor of 1.5 for the stellar masses, ±0.3 mag for the extinctions, 

and factors of ~ 2-3 for the ages as well as for the absolute and s pecific star fo rmation rates. 

" Classification according to the Ha kinematics from SINFONI ( Shapiro et al. 112008; Forster Schreiber et al.ll2006ll2009h . 
' Effective radius, light concentration index, and ratio of minor to major axis (i?e, «, bjd) obtained from single-component Sersic 
model fits to the two-dimensional /fieo band surface brightness distributions of t he galax ies; the best-fit parameters correspond to 
the median of the distribution of results from 500 GALFIT runs as described in iPaper \ and the uncertainties represent the 68% 
confidence intervals around the best-fit value. For BX 528, the two main merging components "SE" and "NW" are resolved in the 
NIC2 data; the results of simultaneous two-component Sersic model fits are = 3.18:0:2!, " = l-20!S:lcj, and b/a = 0.35 _():o3 for 
BX528-SE, and R, = S.STlfj, n = 2.91!!:^, and b/a = 0.51!^:* for BX528-NW. For BX389, the NIC2 data resolve the small 
companion south of the main part of the galaxy; the results of simultaneous two-component Sersic model fits are Re = 5.89!^ 05, n = 
0.32!^-°ot, and b/a = 0.27 ±0.01 for the main body of BX 389, and R^ = 2.7l!g:^^, n = 3.03!?:75, and b/a = Q.IO^oaI for BX 389-South. 
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TABLE 2 

Properties of clumps identified in the NIC2 images 
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" Here, "clum p" ref ers generically to > 3(t peaks on chai'acteristic spatial scales of 1 — 1.5 times the PSF FWHM above the local background in the //[go 
band images (§ 13. U . Clumps #1 and #2 of the major merger BX 528 are part of the SE component while clump #4 is associated with the NW component. 
In BX 663, clump #2 is close to the center of the galaxy and may be associated with a stellar bulge or the AGN present in this object. Clump #4 of BX 389 
corresponds to the small companion to the south of the main galaxy. 

Projected distance (dproj) relative to the geometiic center of the galaxies, and deprojected distance (d) assuming a circular configuration with inclination 
angle corresponding to the best-fit axis ratio from the NIC2 morphology (Table [T). No projected distance is calculated for the southern companion in 
BX 389 (clump #4), which appears to be located below the galactic plane. 

Radius of the circular aperture used to measure the clump fluxes, and aperture correction based on the curve-of-growth from the empirical PSF light 
profile, corresponding to the flux ratio in circular apertures of radius 1" and fphot- 

^ FWHM is the observed direct full-width at half-maximum from the clump radial hght profiles after subtraction of the local background. FWHM" is the 
intrinsic size derived by subtracting in quadrature the observed FWHM and the PSF FWHM. For unresolved clumps, the PSF FWHM is taken as upper 
limit on the size. 

° Estimates of the clump fractional contribution to the total //i60 band emission of the galaxies. is the contribution assuming the clumps are point 
sources, with the local background from the host galaxy subtracted. The J^^^y, fraction is derived from the flux in the aperture of radius rphot. without 
subtraction of the host galaxy background and without aperture correction. These estimates represent lower and upper limits, respectively, as argued in 
§ 13.11 and are listed for comparison with the adopted background-subtracted, aperture-corrected fractions simply denoted J^'^K 
' Background- and aperture-corrected //i6o band magnitudes of the clumps (uncorrected for extinction, in the AB system). 

8 Background- and aperture-corrected equivalent rest-frame g-band luminosity of the clumps (uncorrected for extinction) relative to the Sun's g-band 
luminosity density (which has an absolute Mg,Q = 5.07 mag). 

Stellar mass estimates for the clumps, inferred from their light contributions (background- and aperture-corrected) and the total stellar mass derived from 
the SED modeling of each galaxy. 
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TABLE 3 

H|f,0- AND /gi4-BAND MAGNITUDES AND COLOR-BASED PROPERTIES OF CLUMPS IDENTIFIED IN THE NIC2 AND ACS IMAGES OF MD 41 
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(log[MQL^ g]) 


(10* Mq) 


(mag) 


'814 

(mag) 


log(M,/Lf'-=') 
(log[MoL^ g]) 


(10* Mo) 


1 

2 
3 
4 

ACS-1 
ACS-2 
ACS-3 


0.56 
1.65 
2.96 
4.77 
6.95 
1.91 
6.89 


1.1 ±0.2 
2.3 ±0.2 

3.2 ±0.1 

5.8 ±0.4 
7.0 ±0.1 

2.9 ±0.4 

7.3 ±0.2 


0.10 
0.10 
0.10 
0.10 
0.10 
0.10 
0.15 


26.20 ±0.04 
26.03 ± 0.04 
26.16±0.04 
26.26 ±0.05 
27.16±0.11 
26.73 ± 0.07 
26.46 ± 0.09 


27.45 ±0.02 
27. 16 ±0.02 
27.03 ±0.02 
26.68 ±0.01 
27.62 ±0.03 
27.37 ±0.02 
26.21 ±0.01 


-0 71+°"* 

"•"-0.17 

-0 75+°" 
-0.89±0.17 
-1.15±0.17 

'•'■^-0.20 

-1.02±0.19 
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6.4!^; 
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A 4+-' 

'■"-0.4 

1 9+" 
'■'-0.7 

0.64_o22 


27.08 ±0.21 
26.82 ±0.17 
27.45 ±0.30 
26.59 ±0.13 

> 27.67 

> 27.67 

> 27.49 


> 29.65 
28.23 ±0.10 
28.84±0.17 

27.77 ±0.06 
28.64 ±0.14 
28.53 ±0.13 

26.78 ±0.03 
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-0.65_o23 
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"■ '■^-0.20 
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Note. — i^ieo and I814 photometry and color-based properties of clumps in MD41, where the set of four clumps identified in the NIC2 image (#1 to #4) is augmented with 
three clumps identified in the ACS map but not in the NIC2 map (see Figure[4). The ratio of stellar mass to rest-frame g-band luminosity (uncorrected for extinction) of the clumps 
is computed from the relationship between observed /gu ~^^160 color and Mi,/!]^^^ derived in Paper 1 (see also ^ |4. 1| . The clump stellar masses listed here are obtained from the 
color-based /L^*"^' ratios infeired for individual clumps instead of being based on a uniform ratio across the galaxy as implicitely assumed for the estimates reported in Table [2] 



^ Projected distance (c/pmj) relative to the geometric center of the galaxies, and deprojected distance {d) assuming the clumps lie in a disk of inclination with respect to the sky plane 
given by the morphological axis ratio (TablelTl. 

Radius of the circular aperture used to measure the clump fluxes. 

Properties based on the fluxes measured in the circular aperture of radius rphoi, without subtraction of the local background from the galaxy and without aperture correction. The 
1 a photometric uncertainties are computed from the noise properties of the images. The uncertainties for the M^/U^^^ and results account for those of the magnitudes and the 
accuracy of the k\4—H\(,o versus M^/V^''^ relationship 0.15 dex; see § |4.1.3t . 

^ Properties based on the fluxes measured in the circular aperture of radius rphoi with local background subtraction and aperture conection. For rphoi = 10 and 0" 15, the aperture 
correction is a factor of 2.09 and 1.63, respectively. The 1 a photometric uncertainties reflect the noise properties of the images and account for both the measurement of clump and 
background fluxes; 3 a limits are given for clumps with background-subtracted fluxes in one or the other band that are smaller than the 3 cr photometric uncertainty. The uncertainties 
for the M^/U^^^ and results are obtained in the same way as for those based on the raw photometry. 



TABLE 4 

Hieo-BAND MAGNITUDES, Ha FLUXES, AND STELLAR AGES OF CLUMPS IDENTIFIED IN THE NIC2 AND SINFONI MAPS OF BX 482 

Clump Jpiuj ^ J'' rphoi ^ From ''raw" photometry ^ From background-subtracted photometry ^ 





(kpc) 


(kpc) 


(arcsec) 


(mag) 


(10"" ergs"' cm" 


1 


4.74 


4.8 ±0.2 


0.275 


23.97 ±0.02 


3.92 ±0.05 


2 


3.02 


5.6±0.9 


0.15 


26.02 ±0.05 


0.46 ± 0.03 


3 


4.36 


6.7± 1.3 


0.175 


25.47 ±0.04 


1.14 ±0.03 


4 


1.97 


2.9 ±0.6 


0.20 


25. 14 ±0.03 


0.78 ± 0.04 


5 


3.68 


7.6± 1.0 


0.175 


25.68 ±0.04 


1.09 ±0.03 


Hce-l 


3.69 


5.8± 1.1 


0.15 


26.03 ± 0.05 


0.57 ±0.03 


Hq-2 


6.31 


6.6±0.8 


0.15 


26.41 ±0.07 


0.42 ± 0.03 



F»l Age H,'l„ F^l Age 

(Myr) (mag) (10"" ergs"' cm"^) (Myr) 



55 ±3 24.50 ±0.03 2.41 ±0.05 55!^ 

115!2^ 27.30 ±0.30 <0.15 > 94 

37!^ 27. 19 ±0.29 0.17 ±0.03 88!.J?° 

239!35 26.48±0.16 <0.17 > 570 

25!* 26.39 ±0.14 0.43 ±0.03 54!-]° 

60!|5 > 27.49 0.25 ± 0.03 < 15 

56*\n > 27.49 O.I9±0.03 < 30 



NOTE. — /fi60-band and Ha photometry and stellar age of clumps in BX482, where the set of five clumps identified in the NIC2 image (#1 
to #5) is augmented with two clumps identified in the SINFONI Ha map but not in the N1C2 map (see Figure [6l. The stellar age is inferred from 
the rest-frame Ha equivalent width using Bruzual & Chariot 1 2003) evolutionary synthesis models for a constant star formation rate, as described in 

§EI1 

^ Projected distance (c/pioj) relative to the geometric center of the galaxies, and deprojected distance (d) assuming the clumps lie in a disk of 
inclination with respect to the sky plane given by the moiphological axis ratio (Table |T1. 
^ Radius of the circular aperture used to measure the clump fluxes. 

Properties based on the fluxes measured in the circular aperture of radius rphoi, without subtraction of the local background from the galaxy and 
without aperture coiTcction. The 1 a measurements uncertainties are computed from the noise properties of the images. The uncertainties for the 
ages account for those of the Ha equivalent width measurements, propagated based on the CSF model curves. 

Properties based on the fluxes measured in the circular aperture of radius / phoi with local background subtraction and aperture conection. For rphot = 
0.'l5. 0f'l75, 0. 20, and 0"275, the aperture correction is a factor of 1.92, 1.68, 1.53, and 1.28, respectively. The 1 a photometric uncertainties 
reflect the noise properties of the images and account for both the measurement of clump and background fluxes; 3 a limits are given for clumps 
with background-subtracted fluxes in one or the other map that are smaller than the 3 cr photometric uncertainty. The uncertainties for the ages are 
obtained in the same way as those based on the raw photometry. 



Properties of stellar clumps in z ~ 2 galaxies 



TABLE 5 

Estimates of the contribution of a central stellar bulge component in the SINS NIC2 sample 

DISKS 



Property ' Case 


BX663 


MD41 


BX389 


BX610 


BX482 


J-"""^' (W,60) . . . 

J^"'-" (M^,) Passive M/L 

jrbuige yi^^ jjjjj ggp^ Av = mag 


0.017 

0.018!!!;™ 
0.039 ±0.001 
u.uy/_o.oo3 


0.035 

0.26_o.o2 
0.64_o.o4 

7 -34+0- 10 


0.034 
0.026 ±0.000 
0.063 ±0.000 
0.19 ±0.00 


0.012 
0.011 ±0.000 
0.026 ±0.001 
0.063_o.oo3 


0.013 

0.050 ±0.001 

0.12 ±0.00 
n oo+o.oo 


a jrbuige (/fj^p) is jj^g estimated maximum contribution to the total Huo band emission of each galaxy of a putative hidden 
bulge component with effective radius = I kpc and Sersic index n = 3 (see § 15. 3t. J^"'^" (M*) is the corresponding 
fraction of the total galaxy stellar mass derived from the rest-UV/optical SED modeling. Uncertainties on the latter reflect 
the ranges of fractional contributions when varying the hidden bulge up to 3 kpc, and n between 2 and 4. The range for 



the fraction of light for n = 2 and 4 are within < 0.5% ot the values for n = 3 for all galaxies. 

Different assumptions were made on th e M/L ratio of the hi dden bulge: a ratio corresponding to the median value of the 
massive quiescent z ~ 2 galaxy sample of Ivan Dokkum et"ari (2008), and that of a maximally-old single SSP whose light 
is unobscured or attenuated by the same amount as inferred from the SED modeling of each galaxy (see § l5.3t . 
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BX528 , I BX663 I MD41 

2 = 2.2683 - I 2 = 2.4332 I 2=2.1704 




Fig. 1 . — HST NICM0S/NIC2 //i6o-band maps of the galaxies, probing the broad-band emission around rest-frame 5000 A. The target name 
and Hof redshift are labeled in the top comer of each panel. The color coding scales linearly with flux density from black to white for the 
minimum to maximum levels displayed (varying for each galaxy). Contours are overplotted starting at flux densities ~ 50% of the maximum in 
each image. The FWHM of the effective PSF is indicated by the filled circle at the bottom left of each comer; the angular resolution is 0" 145, 
or ~ 1.2 kpc at the median z = 2.2 of the sources. The angular scale of the images is shown with the l"-long vertical bar next to each galaxy. 
In all maps, North is up and East is to the left. The compact clumps identified in each galaxy are indicated, with their ID number labeled. 
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Fig. 2. — Properties of clumps identified in the NIC2 images of our galaxies, (a) Him band magnitude versus observed FWHM. Dotted 
lines indicate curves of constant surface brightness as labeled, in units of mag arcsec"". (b) Fraction of the total galaxy H160 band light 
measured for each clump versus observed FWHM. (c) Derived stellar mass of the clumps versus observed FWHM. Different colors are used for 
clumps belonging to different galaxies, following the scheme of the labels in each panel. Clumps associated with the southeast and northwest 
components of the merger system, BX 528, are plotted with triangles (upward and downward pointing, respectively), the one corresponding to 
the southern companion of BX 389 is plotted with a square, and the one identified as the central peak of the AGN galaxy BX 663 is marked as 
starred circle. In each panel, the histograms show the distributions of the measurements projected onto the horizontal and vertical axes, and the 
hatched bars correspond to the respective median values (including the upper limits on the observed FWHMs). The vertical dashed line shows 
the FWHM of the PSF. The observed FWHM of the clumps corresponds to the direct FWHM size; for unresolved clumps, the PSF FWHM 
is taken as an upper limit on the observed size. Magnitudes and fractional light contributions are measured in the circular apertures listed in 
T able [2I with the local background from the host system subtracted and an aperture correction based on the PSF profile applied, as described in 
§ 13.11 The stellar mass estimates are computed from the fractional light contributions assuming a constant observed M* /Lhi60 ratio across each 
galaxy. 




0.0 0.5 1.0 1.5 
Intrinsic FWHM"' [kpc] 



Fig. 3 . — Stellar mass versus intrinsic FWHM size of clumps identified in our NIC2 sample. Data are plotted only for the resolved clumps, i.e. 
with observed FWHM > PSF FWHM. Symbols are the same as used in Figure[2] and the histograms and hatched bars show the distributions 
and median values for this set of resolved clumps. 
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Fig. 4. — Comparison of rest-frame optical and UV morphologies for MD41. North is up and east to the left in both images, and a bar 
of l" in length indicates the scale in the images. (Left) NIC2 //leo-band map, with the FWHM = 0"l45 of the PSF indicated by the filled 
circle at the bottom left comer. (Middle) Same as the left panel but for the ACS !8i4-band map, at the original resolution of the data with PSF 
FWHM = 0"099. (Right) Color-composite after accurate PSF matching and registration of the NIC2 and ACS maps; in this RGB image, the red 
channel is assigned to the //i6o band, and the green and blue channels to the igw band. The colors for all images plotted correspond to a linear 
flux density scale. The white cross shows the location of the geometric center of MD 41 as determined from the NIC2 map. Overall, the images 
are strikingly similar, although on small scales there are noticeable differences. In particular, the two brightest clumps in the ACS image on 
the southwest edge of the galaxy are much fainter in the NIC2 image and, consequently, have bluer colors than the bulk of the source. Clumps 
identified in the NIC2 //i6o band image are labeled in the left panel, and those identified in the ACS igw band map that have no counterpart in 
Hido band emission are labeled in the middle panel. 




Fig. 5. — Observed near-IR H\(,o band and optical (8i4 band magnitudes from the HST ACS and NIC2 imaging of MD41. (a) Data for all 
individual pixels across MD41 with a S/N > 3 in both H\(,o and (8i4 bands. Error bars correspond to 1 ct uncertainties on the measurements. 
Pixels that fall within the photometric apertures of different clumps are color-coded according to the scheme given by the labels. Black data 
points show pixels in the interclump regions. Diagonal dashed lines correspond to constant (8i4-//i6o colors, as labeled in the plot; the thick 
yellow-orange line shows the median color of pixels in the interclump regions. The arrows indicate the effects of Ay = 0.5 mag of extinction on 
the magnitudes, and of an emission line contribution in the Hi^o bandpass of 17% derived from the integrated [O III] AA4959, 5007 and H/3 line 
fluxes for MD 41 measured with SINFON I. (b) S ame as panel (a) but for the clump magnitudes. Results based on the "background-subtracted" 
and "raw" photometiy (as explained in § 14. 1. It are plotted as circles and squares, respectively. Error bars (mostly smaller than the symbol 
size) correspond to the formal photometric uncertainties. For clumps undetected in one or the other band, 3 a upper limits on their background- 
subtracted photometry are used. For comparison, the integrated magnitudes of MD41 and those obtained by summing up the fluxes of all 
individual pixels with S/N > 3 are plotted as grey- and white-filled diamonds, respectively. The variations in colors across MD 41 and among 
the clumps are fairly modest for the most part; NIC2 clump #1 and clump ACS-#3 exhibit the most significant differences. The NIC2-identified 
clumps tend to be redder than those identified only in the ACS image. 



Properties of stellar clumps in z ~ 2 galaxies 



25 





Fig. 6. — Comparison of Ha and rest-frame optical morphologies for BX482. North is up and east to the left in both images, and a bar of 1" 
in length indicates the scale in the images. (Left) NIC2 H\bo band image at the original angular resolution of the data, with the FWHM = 0." 145 
of the PSF indicated by the filled circle at the bottom left corner. (Middle) Same as the left panel but for the SINFONI Hq map, with PSF 
FWHM = 0."17. (Right) Color-composite after accurate PSF matching and registration of the NIC2 and Hq maps; in this RGB image, the 
red channel is assigned to the //i6o band and the green and blue channels to Ha. The colors for all images plotted correspond to a linear flux 
density scale. The white cross shows the location of the geometric center of BX 482 as determined from the NIC2 map. The H\fso band and Ha 
emission follow each other closely although on small spatial scales there are some noticeable differences. In particular, the clumpy structure on 
the northwestern side of the galaxy center has less correspondence between the two maps than on the southeastern side. Clumps identified in 
the NIC2 H[(,o band image are labeled in the left panel, and those identified in the SINFONI Ha line map that have no counterpart in //leo band 
emission are labeled in the middle panel. 




I . . . . ■ I J ■y'y'- i . i . i 

0.01 0.1 1.0 10.0 

/hi6o [Wy] /hibo [^Jy] 

Fig. 7.— Observed near-IR H\6o flux density and Ha line flux from the HST NIC2 imaging and SINFONI+AO observations of BX482. (a) 
Data for all individual pixels across BX482 with a S/N > 3 in both /figo band and Ha emission. Error bars show the 1 cr uncertainties on 
the measurements. Pixels that fall within the photometric apertures of different clumps are color-coded according to the scheme given by the 
labels, and black data points show pixels in the interclump regions. Diagonal dashed lines correspond to constant rest-frame Ha equivalent 
widths, translated to stellar ages based on model predictions for constant star formation (see § 14.2.3 I for details); ages in Gyr are labeled next 
to each line (and, from young to old ages, correspond to lV'°^'(Ha) =415, 265, 170, 115, 80, and 60 A). The thick yellow-orange line shows 
the median age inferred for pixels in the interclump regions. The arrows indicate the effects of Av = 0.5 mag of extinction on the fluxes, and of 
the upper limit for the contribution in the H[(,o bandpass of 12% derived from the integrated fluxes of the [O III] AA4959,5007 emission lines 
and the 3(t limit for H/3 measured in BX482 with SINFONI. (b) Same as panel (a) but for the clump measurements. Results based on the 
"background-subtracted" and "raw" photometry (as explained in § 14.2. It are plotted as circles and squares, respectively. Error bars (mostly 
smaller than the symbol size) correspond to the formal photometric uncertainties. For clumps undetected in one or the other image, 3 u upper 
limits on their background-subtracted photometry are used. The integrated fluxes for BX482 and those summed up over all individual pixels 
with S/N > 3 are plotted as grey- and white-filled diamonds, respectively. The model curves are for the same ages as in panel (a). The clumps 
span a comparable range in W"^'*(Ha), hence in inferred stellar ages, as the interclump regions. The NIC2-identified clumps tend to have lower 
ratios, suggesting older ages, than those identified only in the Ha map. 
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Fig. 8. — Properties as a function of deprojected galactocentric distance of cl umps in BX482. (a) Rest-frame Ha equivalent width, computed 
from tiie SINFONI Ha map and tire NIC2 Hyaa band image as described in § 14.2.21 Results based on the "background-subtracted" and "raw" 
photometry are plotted as circles and squares, respectively. Different colors indicate different clumps, as labeled in the plot. Error bars 
correspond to the formal photometric uncertainties. For clumps undetected in one or the other image, 3ct upper limits on their background- 
subtracted photometry are used. The thick yellow-orange line shows the median ^"^'"(Ha) measured for pixels in the interclump regions with 
a S/N > 3 in each map. {b) Same as panel (a) but for the s tellar a ges derived from comparison of the lV"^'*(Ha) measurements to model 
predictions for a constant star formation rate, as described in § 14.2.31 The uncertainties of the ages account for the photometric measurements 
uncertainties. The clump closest to the center of BX482, NIC2-#4, is clearly distinct from the ensemble of the other clumps farther out in 
radius, with its lV"^^'(Ha) implying the oldest age of all clumps. All other clumps have comparable or higher lV"^^'(Ha), hence younger inferred 
ages, than the median for the interclump regions. For different star formation histories (SFHs), the absolute ages can change significantly but 
the relative ages between clumps and interclump regions are qualitatively unchanged as long as the different regions have similar SFHs. 
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Fig. 9. — Properties as a function of deprojected galactocentric distance of clumps in MD41. (a) Observed !8i4-//i60 colors, computed from 
the ACS i8i4 band and N1C2 Him band maps as described in § 14.1.21 Results based on the "background-subtracted" and "raw" photometry 
are plotted as circles and squares, respectively. Different colors indicate different clumps, as labeled in the plot. Error bars correspond to the 
formal photometric uncertainties. For clumps undetected in one or the other band, 3 a upper limits on their background-subtracted photometry 
are used. The thick yellow-orange line shows the median color measured for pixels in the interclump regions with a S /N > 3 in each band. 
(b) Same as panel (a) but for the ratio of stellar mass to dust-attenu ated re st-frame ^-band luminosity derived using the relationship between 
observed /gw -Wi60 colors and M*/L^'^'* at z = 2.2, as described in § 14.2.31 (see also Paper I). The uncertainties of the M/L ratios account for 
those of the photometry as well as for the estimated accuracy of the color-M/L relationship. Clumps closer to the center of MD41 tend to 
have redder colors, hence higher observed M^/Vg"' ratios, which could reflect increasing stellar ages and/or dust extinction of clumps at smaller 
radii. 
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Fig. 10. — Images and light profiles extracted from the NIC2 images of the five disks among our SINS NIC2 sample. Each row corresponds 
to one galaxy, and the panels are as follows. Left: Hibo band maps on a linear greyscale intensity scale. An ellipse is drawn, centered at 
the geometric center of each galaxy, with semi-major axis length, axis ratio, and P.A. corresponding to the best-fit parameters of the single- 
component Sersic model fits of Paper I. The rectangle shows the pseudo-slit used to extract the major axis light profiles. North is up, east is 
to the left. Middle: Radial light profiles extracted in annuli of increasing radius, with center, axis ratio, and P.A. given by the best-fit single- 
component Sersic model ("disk") to the Hieo band 2D surface brightness distribution of each galaxy. The profiles are shown for the data (black 
solid line and filled dots), the best-fit "disk" component (thick orange line), the "maxim al" b ulge component with i?c = 1 kpc, n = 3, and b/a= I 
allowed by the 1 a measurements uncertainties at r < 2 kpc (magenta solid line; see § 15. 3t , and the sum of the "disk" and "maximum bulge" 
(blue solid line). The circularly-symmetric light profile of the PSF is plotted for reference (green dashed line). Right: IVIajor axis light Hieo 
band light profiles extracted extracted along the pseudo-slit of width 0"3 shown in the left panels; negative (positive) offsets relative to the 
center correspond to the south (north) side of the galaxy. Different colors for the lines are the same as for the middle panels. 
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Fig. 11. — Comparison of stellar ages of clumps in BX482 derived from their Ha equivalent width for different assumptions on the star 
foimation history (SFH), as described in Appendix |A] In all panels, the horizontal axis corresponds to the age derived for a constant star 
foimation rate (CSF). The vertical axis corresponds to the age derived for alternative SFHs as follows. Top row: exponentially-declining star 
formation rates with e-folding timescales r = 300, 100, and 30 Myr (labeled "r300," "rlOO," and "r30," respectively, in each panel). Middle 
row: exponentially-increasing star formation rates with the same e-folding timescales (labeled "ir300," "irlOO," and "ir30," respectively). 
Bottom row: "delayed" star formation rates with the same e-folding timescales (labeled "dr300," "drlOO," and "dr30," respectively). Results 
based on the background-subtracted and raw clump photometry are plotted as circles and squares, respectively. Different colors indicate 
different clumps, as labeled in the top left panel. Uncertainties of the ages are computed by propagating the measurement uncertainties on the 
model curves. The thick yellow-orange lines in each panel show the stellar age derived from the median Ha equivalent width of pixels in the 
interclump regions with a S /N > 3 in /f leo band and Ha emission. 



